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Abstract 



Some considerations about the importance of coherence effects for bremsstrahlung 
processes in non-equihbrium dense matter (Landau - Pomeranchuk - Migdal - ef- 
fect) are presented. They are of particular relevance for the application to photon - 
and di-lepton production from high energy nuclear collisions, to gluon radiation in 
QCD transport, or parton kinetics and to neutrino and axion radiation from super- 
nova explosion and from hot neutron stars. The soft behavior of the bremsstrahlung 
from a source described by classical transport models is discussed and pocket cor- 
rection formulas for the in-matter radiation cross sections are suggested in terms 
of standard transport coefficients. The radiation rates are also discussed within a 
non-equilibrium quantum field theory (Schwinger - Kadanoff - Baym - Keldysh) 
formulation. A classification of diagrams and corresponding resummation in physi- 
cally meaningful terms is proposed, which considers the finite damping width of all 
source particles in matter. This way each diagram in this expansion is already free 
from the infra-red divergences. Both, the correct quasi-particle and quasi-classical 
limits are recovered from this subset of graphs. Explicit results are given for dense 
matter in thermal equilibrium. The diagrammatic description may suggest a formu- 
lation of a transport theory that includes the propagation of off-shell particles in 
non-equilibrium dense matter. 



permanent address: Moscow Institute for Physics and Engineering, Russia, 115409 
Moscow, Kashirskoe shosse 31 

Electronic mail: J.Knoll@gsi.de, voskre@rzri6f.gsi.de 



GSI-Preprint 95-63, subm. to Ann. Phys., (hep-ph/95****) 



1 Introduction 



The importance of colierence time effects on the production and absorption of field 
quanta from the motion of source particles in non-equilibrium dense matter has first 
been discussed by Landau, Pomeranchuk, Migdal (LPM) [1,2] (and many others 
later) in the context of bremsstrahlung from ultra-relativistic electrons undergoing 
multiple rescatterings on Coulomb centers. The first successful measurements of the 
corresponding suppression of bremsstrahlung have been carried out at the Stanford 
Linear Accelerator Center very recently [3] . With this paperj^ we like to supplement 
some quite intuitive and also formal considerations, which illustrate the nature of 
production and absorption processes in a dense matter environment. The subject is 
of quite general nature and applies to many physical problems, where either a source 
couples weakly to a wave field or for the proper determination of local gain and loss 
terms in quantum transport. Examples are the application to photon, or di-lepton 
production from a piece of dense nuclear matter or hadron gas formed in high energy 
nuclear collisions, for gluon or parton radiation and absorption in QCD transport and 
its practical implementation in parton kinetic models (such problems are discussed, 
e.g. in [5,6]), to neutrino and axion radiation from supernovas and neutron-star 
matter (see [7,8,13,10]), for the soft phenomena in quantum cosmological gravity 
(see [11]) and also for many condensed matter phenomena, as particle transport in 
metals and semiconductors, radiation in plasma etc. (see [12,13]). To be specific, 
however, we take the example of electrodynamics, considering photon production 
from a piece of nuclear matter, but when appropriate comment on other cases. 
Since throughout the paper we discuss the corresponding proper self energy of the 
produced particle, all considerations given also apply to gain and loss terms of other 
particles in non-equilibrium dynamics. 

In the context of high-energy nucleus-nucleus collisions [6] for example, it became 
quite apparent over the last years, that a justification of QCD transport (e.g. in terms 
of a parton kinetic picture) calls for a proper understanding of all soft processes. 
Well known is the Rutherford singularity in scattering cross-sections of interactions 
mediated by the exchange of zero mass quanta (photons/ gluons). In dense matter 
the exchanged quantum acquires a finite real mass due to Debey-screening. Singu- 
larities are also encountered in absorption or radiation processes (bremsstrahlung). 
Induced by free scattering the rates diverge at vanishing four-momentum q of the 
radiated quantum, due to the infinite time scales used in the quasi-free approxima- 
tion. In dense matter, however, due to the finite free propagation time t^ow between 
successive collisions these rates become regular. In particle physics context, most 
of the papers on the LPM effect discuss the bremsstrahlung of some fast charged 
particle, such as hadron, quark or gluon, which traverses a dense hadron gas or 
quark-gluon plasma, e.g. see refs. [14,5], where the role of the matter is reduced to 
infinitely massive scatterers. In reality all the particles in dense matter which couple 
to the radiated field should be treated on equal footing. Effects of the finite mean 
free propagation time on photon and gluon radiation have been considered e.g. in 

^ a brief report of these results is given in [4] 
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refs. [14,16-18,5]. 



While the problem can be quite simply and intuitively formulated and solved in the 
classical limit, where a classical source couples to a wave field, e.g. classical charge 
particles couple to a Maxwell field, considerable conceptual difficulties arise for the 
very same problem, if the source is described as a quantum many-body system. 
In fact common standard techniques, like perturbation theory or the quasi-particle 
approximation (QPA) have serious limitations to describe the production and/or 
absorption rates over the whole range of energies and momenta, as they completely 
fail in the soft limit. 

Starting from a quantum many-body formulation in terms of Green's functions most 
derivations of transport descriptions employ two essential approximation steps: i) a 
gradient expansion and ii) the QPA. For simplicity we concentrate on the defects 
of the QPA in this paper. In the QPA, which is a consistent approximation scheme 
for low temperature Fermi liquids (Landau - Migdal, see [19,20]), all particles in 
the medium are treated on-shell with a well determined energy-momentum relation 
(dispersion relation) which follows from the real part of the retarded self energy 
of the particle. To be specific in this notion, we use the term "on-shcH", when the 
particle follows a sharp energy-momentum relation. Thus the quasi-particle poles of 
the retarded Green's function lie just infinitesimally below the causality cut along 
the real axis in energy. The corresponding approximation scheme in terms of these 
on-shell states, which have infinite life time, is called quasi-particle approximation 
(QPA). Due to interactions in dense matter the damping of the quasi-particles may 
become important, the corresponding "quasi-particle" poles of the retarded prop- 
agators move into the unphysical sheet below the real axis. As a consequence the 
mass spectrum of the particles is no longer a sharp delta function but rather ac- 
quires a width r, and one talks about "off-shell" propagation. In that case one has 
to leave the standard description in terms of stable single particle states and em- 
ploy quantum propagators (Green's functions) with continuous mass distributions. 
Landsmann [21] has coined the notion of "non-shell particles" in this connection. 
One thus comes to a picture which unifies resonances which have already a width 
in vacuum due to decay modes with the " states" of particles in dense matter, which 
obtain a width due to collisions (coUisional broadening). 

The theoretical concepts for a proper many body description in terms of a real time 
non equilibrium field theory have already been devised by Schwinger, Kadanoff, 
Baym and Keldysh [22-24] in the early 60'^*. First investigations of the quantum 
effects on the Boltzmann collision term were given by Danielewicz [25], the princi- 
ple conceptual problems on the level of quantum field theory were investigated by 
Landsmann [21], while applications which seriously include the finite width of the 
particles in transport descriptions were carried out only in recent times, e.g. [12,25- 
35]. For resonances, e.g. the delta resonance in nuclear matter, it was natural to 
consider broad mass distributions and ad hoc recipes have been invented to include 
this in transport simulation models. However, many of these recipes are not correct 
as they violate some basic principle like detailed balance [32] , and the description of 
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resonances in dense matter has to be improved. The present study also gives some 
hints on how to generahze the transport picture towards the inclusion of off-shell 
propagations in dense matter [27,30,31]. 

In this paper we illustrate the practical implications of such non-equilibrium concepts 
at the example of particle production from the dense matter dynamics. Thus all 
source particles never reach an asymptotic state and naturally have a continuous 
mass spectrum. In sect. 2 we derive the basic formulas for the rate of bremsstrahlung. 
The classical and general quantum mechanical expressions, the latter in terms of 
non-equilibrium Green's functions and self energies, are derived for the case of non- 
equilibrium dynamics. In sect. 3 we concentrate on the description of radiation from 
classical sources. We start with the bremsstrahlung from a classical diffusion process, 
and subsequently derive the photon spectrum for a classical random walk (Langevin) 
process in terms of a completely regular multiple collision expansion. The low energy 
behavior is discussed and pocket correction formulas for the in-matter radiation 
cross sections are suggested in terms of standard transport coefficients. Also finite 
size corrections are obtained. Then in sect. 4 we use the non-equilibrium Green's 
function formalism, see [22-25,36,29], and formulate diagrammatic resummations 
where all quantities are expressed through physically meaningful terms. We show 
how infra-red convergent results can be obtained through the account of the finite 
damping width and discuss the QPA (sect. 4.5) and quasi-classical (QC) (sect. 5) 
limits from the corresponding infinite series of diagrams. In sect. 6 the lowest order 
loop diagrams for the production rate from a piece of equilibrium dense matter are 
analyzed in the quantum case, both at high and low temperatures. Conclusions and 
perspectives are given in sect. 7. Some formal details are deferred to the Appendix. 

We use rational units h = c = 1. Whenever the behavior of some quantity is discussed 
in the classical limit {h ^ 0), h will be given explicitly. 



2 Basic Formulas for the Rate of Bremsstrahlung 

If the source system couples only perturbatively (to lowest order in e^) to the elec- 
tromagnetic field, the production or absorption rate of photons can be formulated 
using standard text book concepts in terms of Fermi's golden rule. The correspond- 
ing transition amplitude is given by the electromagnetic current operator between 
the initial and final states of the source. For dense matter problems it is more ad- 
vantage to use a more general concept, where the local production and absorption 
rates are expressed through the current-current correlation function{^ 

^ For the description of coordinates and momenta we use the following conventions: num- 
bers 1, 2, etc. abbreviate space-time points xi = (ti,xi), etc.; for two-point functions 
coordinate means are x = (xi + X2)/2 = (t, x), relative coordinates: = xi — X2 = (t, 0; 
the corresponding four vector Wigner coordinates are {x;q) = {t,x;u;,q) for the photon 
and (x; k) = (t, x; e, k) for the particles of the source. Whenever advantage or necessary we 
shall swap from one to another or even to some mixed representation, just changing the cor- 
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(/^(2)i''(l)) for production, and (/(2)/t(l)) for absorption. (1) 



Although j{x) = j{x,t) is a hermitian operator we distinguish between j and 
in order to designate the photon creation and annihilation vertex, respectively. The 
bracket (...) denotes a quantum ensemble average over the source; quantum states 
and operators are taken in the Heisenberg picture. 

With reference to the description of non-equilibrium systems, where it is advantage 
to use real-time non-equilibrium field theory concepts, such as the Schwinger - 
Kadanoff - Baym - Keldysh technique [22-25] we introduce the following notions 

47r(/^(2)/(l)) = -in'^'^-+(l;2), 47r (/(2)j'^t(i)^ = -in^'^+-(2; 1) , (2) 

which relate the correlation functions to the proper self energies n~+ and n+~ of 
the photon, which are responsible for gain and loss (c.f. sect. 4.3). Throughout this 
paper we use the {— ,+} notation, defined in detail in sect. 4 in the convention of 
ref. [36], chapt. 10. 

In this formulation the production term for the phase space occupation nj{x,q,t) 
(Wigner density) of on-shell photons per space-time volume d'^x = d'^x dt, and per 
energy-momentum cu — q volume, d^g = dcud^q, with polarization rj — {77^} is given 

by 

d\{x, q, t) = -i77^77,n'^'^-+(x; q) (1 + n^{x, q, t)) - a;|)d^xd^g , (3) 
where cu^ is the photon on-shell energy, and 

- iU^'^'+ix; g) =47r I d^^e^^^ (fH^ - + m) (4) 

denotes the space-time Wigner transformation of the auto correlation function (2). 
This local gain term is the on-shell version of a general quantum transport con- 
cept (Kadanoff-Baym equation [23], c.f. (46), sect. 4). It hkewise applies for virtual 
photons (e.g. dilepton production), replacing the on-shell 5-function in (3) by the 
corresponding off-shell photon spectral function. 

The above expressions are the space and time-dependent version of the more familiar 
golden rule for quantum transitions between exact stationary eigenstates. Note that 
by definition, c.f. (4), —irji^rj^U.'^^ is a real quantity; if integrated over phase-space 

responding arguments of the functions; for space -, or space-time independent systems we 
drop the argument x or x, respectively; e.g.: = 11(1; 2) Il{x; q) — n(a;, q) ^ n(r, q); 

the latter two in space-time homogeneous systems. For simpUcity the polarization indices 
fi and or i and k for the spatial part of the tensor structure of 11'^'^ will not always be 
given in later equations. 
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volumes AxAq large compared to h it becomes positive and serves as a production 
rate. 



Such a formalism has been applied in many cases employing the QPA for the equi- 
librium Green's functions, c.f. refs. [39,10]. However, the general formalism allows 
to go beyond this limit and to account for the finite damping width of the source 
particles due to their finite mean free path, which is the main topic of this paper. 

Therefore the current-current correlation function is the central quantity of interest. 
In graphical form it is determined by the proper self energy diagram of the photon 




in-+ = --{-in}-- (5) 



which sums all one-photon irreducible self energy diagrams[^ The dashed lines relate 
to the photon, while the interior area (—ill) symbolically denotes the exact inclusion 
of all strong interactions among the source particles. 



2.1 Analytical properties and constraints 

The self energies for gain and loss obey some analytical relations that follow right 
from the definitions (2,4), like 

-in^^-+(l;2) = (-in'^^~+(2;l))* 

which implies that —i7]fj_riiyll'^'^{x; q) is real. Production and absorption parts obey 

W^~+{x]q) = W^+'{x]-q). 

Integration over uj projects onto equal time properties. Of particular help for the 
discussion of soft processes [16,17] are the following energy weighted dipole (g = 0) 
sum-rules (e.g. [37]) 



-i / ^uj''-^U-+{u,q = 0-t,x)d^x = Sn= < 



4:-K(^j\t)J^{t)) for n = 2 

-27ri([L)*(t), J'=(t)]) for n = l (6) 
Att {D\t)D^{t)) for n = 



which are valid also in the general non-equilibrium case. Here i,k ^ {1,2,3} de- 
note the spatial components. The r.h. expressions are given by the space integrated 

^ To order naturally all diagrams are one-photon irreducible; however for the applica- 
tion to the production of particles with a larger coupling constant, e.g. for gluons with 
coupling constant g, also diagrams to higher order in g are relevant and one then has to 
discard diagrams which are one gluon line reducible. 
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currents and dipole moments 

r{t) = J d'xf{t,xy, D\t) = j d?xx'f{t,x). (7) 

While the n — 2 sum-rule directly follows from definition (4) and applies to any 
current even non conserving, the other two use current conservation and partial 
integrations, known as Siegert's theorem (c.f. [38]), and therefore also require that 
the system has a finite space extension. (If there are no long range correlations, the 
n = 2 and n = 1 relations can also be used for infinite matter, if properly taken 
per volume). For non-relativistic currents the commutator in the n — 1 (Thomas 
- Reiche - Kuhn) sum-rule just becomes the sum of square charges in the system 
i[D^{t), J*(t)] = Y^e-l, where u labels the constituents. 

For systems in thermal equilibrium production and absorption follow the detailed 
balance relation (Kubo - Martin - Schwinger [15]) 

n-+(g; x) = n+-(g; x)e-'"/^ (8) 

where T is the temperature. They allow to write the l.h.s. of the sum-rules as half- 
sided integrals, e.g. integrating only the production rate 

- i I ^a;"-^n-+(u;, 0; t, x) (l + (-l)"e-/^) d'x = S^. (9) 



These rules have been used to estimate the validity of the quasi-free scattering 
prescription in kinetic models [16,17]. In the classical limit, where hw (c.f. sect. 
3. below) the l.h.s. of the n = 1 and n = 2 sum-rules coincide in equilibrium, apart 
from a factor T/2, and the ensuing identity 25*2 = TSi is a disguised form of the 
classical equal partition theorem. 

A further consistency check for diagrammatic elements which determine the self 
energy can be given in terms of Ward identities in the case of conserved currents. 
Since the space-integral of gives the conserved total charge Z of the system, one 
may also use that the space integrated density-density correlator is constant in time, 
i.e. 

J d^a;id^a;2n°°-+(l; 2) = J d^a;n°^-+(r, q = 0;t,x)= AttZ^ = const, (10) 

which applies even in non-equilibrium cases. For isolated systems the motion of the 
center of mass leads to no radiation. Therefore one normally introduces effective 
charges for the different kind of particles of the source in the standard manner such 
that the total effective charge vanishes Zgff = 0. 
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3 Radiation from Classical Sources 

In this section we discuss two examples wliicti treat tlie source as a classical system 
coupled to a Maxwell field. This limit just amounts to evaluate the current-current 
correlator on the classical levelQ. We discuss the radiation caused by a single charged 
particle (the source), which stochastically moves in neutral dense matter. The motion 
of the source is described (a) by mesoscopic transport (diffusion process) and (6) 
by a microscopic Langevin process. Since these examples represent the QC limits of 
the corresponding quantum field theory cases, we carry on the discussion in terms 
of the photon self energy n""*". 



3.1 Diffusion Process 

The motion of a non-relativistic source particle is assumed to be described by a 
time dependent phase-space distribution f{x,v,t) in space and velocity with con- 
vective current density j{x,t) = ej d^v v f{x,v,t). For standard dissipative media 
in equilibrium the velocity autocorrelation function (integrated over space) decays 
exponentially in time 

(^v\t)v''{0)) = i 5^'=e-^-l"l, (11) 

where is the relaxation rate which is supposed to be approximately constant on 
the relaxation time scale l/T^. It relates to the spatial diffusion coefficient D via 
Einstein's relation 

1 °° 1 
D = - jdr {v{r)m) = ^ {v') ■ (12) 

Compared to the infra-red divergent quasi-free result oc l/u'^ (c.f. eq. (28) below) 
this form of the correlation renders the photon self energy 

- in,l+(c., q = 0) = Wpo (v'v^) = 4vreVo4^^^' (13) 

regular at four momentum g = 0. It is determined by mesoscopic transport proper- 
ties, namely by the diffusion coefficient D and relaxation rate r^,; po is the spatial 
density of the charged particles. 



^ One has to realize that a classical photon carries no energy in the quantum sense, i.e. 
hio ^ and the energy is given by the electromagnetic fields. 
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Both f(x, V, t) and the autocorrelation function can be obtained in closed form, if 
the time evolution of /, and the propagation of fluctuations 6f are governed by a 
standard (non-relativistic) diffusion process (Fokker-Planck equation) 



d_ 
di 



f{x,v,t) 



(14) 



In the equilibrium limit (t oo) the distribution attains a Maxwell-Boltzmann 
form 



feq{x, V) = pofeq{v) = pO {^TtDT^) exp 



2Dr, 



m 



-(e(^)-M)/r 



, (15) 



where T — m {iP) /2> — mDV^ and p are the equilibrium temperature and chemical 
potential and e{v) is the energy of the particle. 

At r = we consider an initial fluctuation 6f{x,v,T = 0) = 6'^{x)6'^{v — vq)- 
Its propagation in the equilibrated matter is also governed by the Fokker Planck 
equation (14). By a Gaussian ansatz for the Fourier transform of this fluctuation 
Sf{q, y,T) — J d^xd^vf{x, v, r) exp[—iqx + iyv\ the time-dependence can be obtained 
in closed form as 



Sf{q, y, t) — exp —A + iBy — Cy , where 
C^^(l- e-^r^-) , B = voe-^^^ - iqD (l - e'^^^)' , 



(16) 



A^iq 



T 

I dT'B(T') 



2r. 



2r^T — e 



+ 4e' 



-3 



iq vo 



(e-^^^ - l) 



This fluctuation Sf is the conditional probability which determines the time-depen- 
dence of the current autocorrelation function. With four vectors {v^} — {l,v} and 
{5^} = {1, B} one can express the full correlation tensor in the mixed r, q repre- 
sentation as 



-in^r''(T,g1 =47r J d'xe-'^^{f(x,r)j^{6,0 



I eq 



47reVoexp j-^ (r,|T| + e"^-!^! - l) 



(17) 



X < 



{v^v")^^ e-'^-l^l - D'^q^'q'' (e"'^-!^' ~ ^) \ ^ ^ ^' 

iq'^D (e-^-\''\ - l) sign(T) for = 0, i/ e {1, 2, 3} 

1 lox jJL — V — {) 
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with A and S as a function of r, q and Vo = ^ from (16). Here the ensemble average 
(• • •)eg equihbrium distribution f^g keeps only even moments of v with 

({7^) = SDTx- The result complies exactly with current conservation, i.e. one verifies 
doW^" + iqkU'"' = 0. 

For transverse photons terms proportional to q^q'^ drop. The corresponding spatial 
part of the tensor is shown in fig. 1, right part. This correlation function decays 
exponentially as ~ e"'""''^ at g*= 0, and its width further decreases with increasing 
momentum q = \q\ due to the increase in spatial resolution. The left part shows 
the corresponding density-density correlation {/i — u — 0), which decays only for 
non-zero momentum, due to charge conservation. 



Density corr. function Current corr. function 
1.2 I 1 1 1 1 1 1 1-2 I 1 1 1 1 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Fig. 1: Density-density and current-current correlation functions, — ill™ ""'"(r, ^) and 
— in^2^ '"(r, g), normalized to the values at r = as a function of time r (in units of 
l/Fj;) for different values of the photon momentum = 3fe^r^/< v'^ > with A; = 0, 1, 2, 3. 

The remaining time Fourier transformation gives the cu, g-dependence of the pho- 
ton self-energy. It can be expressed in terms of the incomplete gamma function. 
Straightforward expansion in powers ] (Dg^/rj;)e~'"^l'^l [ leads to 




in,l+(a;, q) = Aire'po [v'v^)^^ exp [D <i'/T, 

-D^Y 2{n+l)V^ + 2D(f 
n\ \ I ((n+l)r, + Dg^)' + cu2 ^ 

for transverse photons. Since the correlation functions are properly determined from 
the time structure of the source, they comply with the n — 2 (and n — 1) sum- rule 
constraints. 

There are two limiting cases where simpler analytical forms can be obtained: i) 
at small momentum transfers where eq. (18) can be expanded in powers of q and 
rewritten as to provide a propagator type form and ii) for large momentum transfers 
where from the exp{. . .} part in (17) a short time Gaussian behavior emerges. Thus, 
for small momentum transfers one finds 
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lim |-in,l+(u;,g)| = 47reVo^'^^*^) Sfx^ 



3(4r2+a;2) 

which generahzes the relaxation result (11,13) to finite q. On the other hand for 
large momenta one realizes that 



lim 
lim 



iU;+{r, q)\ = 47reVo (^'^')^^ exp [-D ^r,ry2 



and therefore 



-in,l+(u;,g)J=47reVo(^^V'=) 



27r 



UJ 



exp 



69 V -C g^r. 



(20) 



27r|g| 



where obviously the essential contributions come from velocities which satisfy the 
Cherenkov condition \v\ f» a;/|g|. This limit is independent of the relaxation rate r^; 
and coincides with the quantum one-loop diagram result in the corresponding large 
|g| limit, as we shall see in sect. 6. 

Although the above expressions give the exact solution of the mathematical problem 
posed in this section, its physical interpretation has to be done with some care for the 
following reason. The equilibrium source distribution contains velocity components 
that exceed the speed of light. Therefore for the physical result mistakes of the order 
of exp [— 3/(2 (f^))] = exp[— m/(2T)] are expected. This restricts the application to 
non-relativistic sources and for large q to space-like photons, where |g| 3> a;. 

For systems with given fixed mean-square velocity (f ^) = const, the exact classical 
on-shell rate (3) at |g| = (x; evidently scales as a function of cj/F^. It properly vanishes 
at a; = and at infinity. It is important to note that the rate has an upper bound 
of ~ ^e^po (f^^)) and indeed attains its maximum value around cu ~ F^., which is 
collision-rate independent. For simplicity we quote the closed form obtained in the 
non-relativistic limit (19), which coincides with the dipole limit. There 



d^ 



d^xdudt 



4 2 / 2\ ^/r^ 



4 2 / 2 



Tx/oJ for a; ^ F-r 

I fora; = F, (21) 

oj/Vx for UJ <^ Fx- 



One realizes that the ultra-violet part of the spectrum uj ^ Tj. behaves as intuitively 
expected, fig. 2: the rate grows proportional to the relaxation rate, until it saturates 
around uj ^ F^- For the soft part a; -C F^;, however, the rate becomes inversely 
proportional to the collision rate! The higher the collision rate the more suppressed 
the spectrum. In order to illustrate the non-perturbative character of this soft be- 
havior supposes F3. oc g^, where g is the strong coupling constant of the source 
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system. One sees that indeed the low-a; part with U oc e^cu/T^ oc e^uj /g^ represents 
a genuine non-perturbative result in g, while the large a;-part, where H oc e^g"^ /u, is 
well described perturbatively. 



Photon rate (diffusion model) 

r-;—! J TV"' 1 1 1 1 1 1 1 1 1 <■ 

1.4 - ; \ \ 



1.2 - 
1.0 - 
0.8 - 




50 100 150 200 250 300 350 400 



« (MeV) 

Fig. 2: Rate of real photons dN / [Aiodt) in 
units of 47re^ (i;^) /3 for a non-relativistic 
source for =50,100,150 MeV; for com- 
parison the IQP results (dashed hues) are 
also shown. 



Correlation funct (Langevin model) 



r 1 1 1 1 r 


— 1 1 1 1 






\ n=0 




" n=1 \. 




n=2 

7 _...--"'_n=3.V^ 
— 1 r 1 1 


1 i —1 -""l " 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

Fig. 3: Current correlation function for the 
first terms n = 0, 1, 2, 3 of the Langevin re- 
sult, eq. (22) and the total sum (S) for the 
case that 

r./r = {{Vm - Vm+l?) /{2{^)) = 1/3. 



3.2 Microscopic Langevin Process 



For a later comparison with quantum diagrams in sect. 4 we should look into the 
corresponding microscopic picture of classical propagation. There one can consider a 
classical random process (Langevin process), where hard scatterings occur at random 
with a constant mean collision rate F. These scatterings consecutively change the 
velocity of a point charge from Vm to Vm+i to ... (in the following subscripts m, 

n, and I refer to the collision sequence, while superscripts i,k & {1, 2, 3} specify the 
spatial components of vectors and the self-energy tensor). In between scatterings the 
point charge moves freely. For such a multiple collision process some explicit results 
can be given. They all refer to the case of vanishing photon momentum g = (dipole 
approximation) and therefore apply to non-relativistic sources where (v^) -C F^ 
or to dilepton production, for example, since only the time structure is well known 
in this case, while the space structure would require an integration of the random 
classical paths. 

For such a Langevin process the modulus of the autocorrelation function takes a 
simple Poissonian form (fig. 3) 



- in,l+(r, q = 0) =47reVo {v\r)v\Q)) 

CO IP In 

= 4.eVoe-l-lE^(->^-^nL■ (22) 
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Here (• • •)m denotes the average over the discrete coUision sequence {m}. This form, 
which one writes down intuitively, directly includes what one calls damping and 
therefore corresponds to a resummation description in the quantum case. The cor- 
responding perturbation theory result is obtained through an expansion in powers 
of r 

- in,r (r, g = 0) = 47reVo E ^ E(-l)' ; (<^^+^>^ ' (23) 

n=0 1=0 \V 

which for dimensional reasons is also a power series in |r| in this case. If the 
^rn^m+Ti) expectation values are replaced by unity in (22) or (23), one obtains 
the 00-component of 11 which becomes constant in time in line with (10). 

The time Wigner transform of (22) determines the a;-spectrum at vanishing q 

- iU:,^., g = 0) = Wpo E {->t^n)^ Re I '^(jf^r^r^ } • (24) 

This is a genuine multiple collision description for the photon production rate in 
completely regular terms due to the (cu^ + F^)*^ form of all denominators. Each term 
is regular, since right from the beginning one accounts for the damping of the source 
particle because of its finite mean time l/F between coUisions. The result (24) still 
accounts for the coherence of the photon field, now expressed through the correlations 
{vmVm+n) in arising from the sequence of collisions. Note in particular, that, although 
the total expression is positive, the n > terms can be negative since they describe 
the interference of the radiation arising from different propagation segments of the 
source particle. Thus, the terms in (24) define partial rates, which later (sect. 5.2) 
will be associated with specific self energy diagrams. 

As already mentioned, the g*-dependence of the self energy cannot be given in closed 
form in general apart from the n — term (c.f. with n — term from eq. (24)) 

-in^^(..a^4.eVo( (^f;^;p )^. (25) 

It shows the typical Cherenkov enhancement at u! = qv. At this level one may be 
tempted to associate this (n = 0) term with the relaxation ansatz result (13). This 
however is only true if < Vm Vm+n >m= forn 7^ 0, an approximation recently used 
in refs. [14]. In the general case velocity correlations between successive scatterings 
exist, and there will be a difference between the microscopic mean collision rate F and 
the mesoscopic relaxation rate F^.. Still, for systems, where the velocity is degraded 
by a constant fraction a per collision, such that ■ Vm+n)m ~ ("^m ■ '^m)^, one 
can resum the whole series in (24) and thus recover the relaxation result (13) at 
^ = 0. The macroscopic rate is then determined by the microscopic scattering 



13 



properties through r^, = (1 - a)r, or 2 {ivm)^)m = ((^m - Vm+i)^)^ F. This clar- 
ifies that the diff^usion result (18) represents a resummation of the random multiple 
collision result. 

The foUowing relations show different reformulations and hmits of the Langevin 

result (24). For instance the invariance of (24) is not directly visible, since absolute 
velocities enter. Still the perturbation expression (23) can be rewritten, such that 
except for the zero order term, which drops out in the Fourier transform, only 
velocity differences appear 



-inei'^(T, g = 0) = 47reVo i {v'A 



iFrl 



iFr 



n n—2 



n=2 1=0 



. I > 



m+l+1 '^m+l+2 



(26) 



Terms of lowest odd order in \r\ determine the asymptotic large u (ultra violet) 
behavior of the spectrum 



hm [-in 



2n 2n-3 



+ E M E(-i)' 

n=2 /=0 




2 \('^m '^m+l)('^m ^m+l) 



Vlr, t';„+i)(t'^+/+i V. 



k 

m+i+2 



(27) 



Apart from the mean collision time ~ 1/F this is an expansion in powers of (F/a;)^ 
and therefore represents the perturbation expansion result for the classical source 
(F representing the interaction, while 1/uj relates to the intermediate propagator). 
This perturbation expansion (27) is interesting since it already displays the main 
problem: While for a; ^ F the series converges, if higher order correlations cease 
sufficiently fast, there is no hope to ever recover the correct result (24) for a; <^ F. 
This is so, since i) this series is necessarily divergent (it has to recover l/(a;^ + F^) 
by a power series in F), but also it misses the knowledge on the even powers in r in 
eq. (26) which essentially determine the soft behavior. 

The first term in (27) represents the incoherent quasi-free production rate which for 
finite q at given polarization 77 reads 



- "^VnVv^iqp (^, = 47re poF 



m+l 



qkjn qkm+i 



(28) 



It carries the known divergence at the soft point g = 0, c.f. dashed lines in fig. 
2 and 4. In conclusion: the commonly used IQF prescription fails for soft particle 
production. 
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3.3 Finite Size Corrections 



For systems of finite spatial extension and conserved currents one can consider the 
n = sum rule. It demands that U^'^ has to vanish at least quadratically 

with q —>■ 0. This property survives in the classical limit, where formally ^ ^ 
and the spectrum becomes continuous. Thus, the term of zero order in uj given by 
I^oD has to drop and the ensuing low energy part of the spectrum 

(24) starts quadratically in u 



hm [-iU:+{uj,q = 0) 

{fitin + 1)(« + 2) (<.,U„>,„ + O (f) ) ■ '2^' 

The simple relaxation ansatz (11) does not fulfill this finite size condition, since it 
ignores long term anti-correlations on the scale of some recurrence time l/F^ec- That 
is the time, where on the mean the center of charges returns to the same pointQ. 
This defect of the relaxation ansatz may be cured by a more general form which 
includes such an anti-correlation, e.g. 



The extra parameter Tree can be determined such that the spectrum fulfills the n = 
dipole sum-rule (6). For larger systems one infers that Tree ~ 2 (f ^) /{{x^) F^,.), where 
(x^) is the mean-square extension. 

For small systems both time scales become comparable and the self energy tensor 
attains the form 

- in,l+(u;, g = 0) = 47reVo {v'v^) .J^_^"p)2 = (31) 



This form has been found in the one-dimensional model of ref. [17], where due to 
the dominance of back scattering both time scales merge. 



not to be confused with the Poincare recurrence time, which is of no relevance here. 
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3.0 



Correlation function (MonteCarlo) 




Fig. 4: Current correlation function 
— in^-*- ^{uJ,q = 0) for a random colli- 
sion sequence limited to a finite size in 
space such that (x^) ~ 10(f^)/r^ and 
= O.SF: full line from Monte Carlo cal- 
culation of the phase integral (32); dashed- 
dotted line from the analytic Langcvin re- 
sult (24); dotted line includes finite size 
corrections (30); dashed line from quasi- 
free scattering prescription. 
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3.4 Comparison with Monte Carlo Evaluation of Amplitudes 

For illustration wc like to present a simple model result and compare it to a Monte 
Carlo method, where amplitudes are calculated by considering the phase of the 
photon field along the classical orbits. Thus, one evaluates 



along the random straight sections of the classical paths in the cascade model. This 
method also permits to illustrate the finite size corrections discussed above. The 
dotted hne in fig. 4. shows the Monte Carlo phase integral result (32) from an 
oriented random walk, compared to the simple relaxation result (13), dashed line. 
The example has the property that the motion is limited to a finite space with 
(x^) ~ 10 (f ^) /F^ and F^; = 0.8F. The sharp dip at = and the little over shoot 
around uj ~ 0.3F are due to the finite size of the system. The full line gives the 
relaxation result including the finite size corrections (30). 

Evidently the Monte Carlo method is a highly unreliable due to the strong cance- 
lations of terms that are randomly generated. The precision in fig. 4 is obtained 
with 200 cascade runs where each path has about 10'^ collisions (simulation codes 
have by far less statistics!). The analytical result (24) has significant computational 
advantages, provided the random process is of this form. For the same precision 
already a single representative path with about 10^ collisions is sufficient, while for 
the relaxation ansatz one only has to determine the relaxation rate, i.e. a simple 
moment in time. 




(32) 
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3.5 Infra-red Divergences, Current Conservation, Gauge Invariance and Identities 

In the above multiple collision description with damping (22) all terms have a finite 
range in time. Thus, they all are void of infra-red divergences and so is any limited 
sum of terms of them. This is the desirable feature that we are aiming at. In addition 
the example above illustrates that in many cases only a few rescattering terms are 
necessary in order to properly recover the correct result for the transverse part of 
the correlation tensor both at small and large cu. In fact the number Ux of required 
rescattering terms is given by ~ r/Fj.. 

The picture is more subtile for the longitudinal components of the tensor, since some 
particular integrals are conserved (time independent) as they related to the total 
charge, c.f. eq. (10), sect. 2.1. Thus — in°'^~+(g = 0, r), c.f. fig. 1, has to be constant. 
For the Langevin case this identity holds, since 

E ^^-^ - 1- (33) 

n=0 "'■ 

However, for any finite number of such terms this identity can never be recovered 
exactly for all times, since the relaxation time of total charge is infinite! 

Quite often one is interested in a solution of the problem up to a maximal time t^ax, 
and one likes to request charge conservation on the correlator level to be maintained 
only within this time span or within the corresponding ti;-range limited to a; > 
i^min — ^/tmax- With this limitation a finite number of terms ~ ^tmax ~ /^min 
is required for a proper description also of the longitudinal parts of the correlation 
tensor at g = 0. 

The above features have to be contrasted with perturbation theory (c.f. eq. (23)), 
where the correlation function can be expressed by a power scries in F in this case 
and therefore leads to a power series in r. Here the zero order term is finite (and 
trivially constant in time), while the higher order terms cancel out for each given 
order, c.f. (23). Thus, in perturbation theory current conservation is maintained 
order by order also for the classical correlator (as one is used to from the quantum 
case). The price to be payed is that in perturbation theory the infra-red properties 
are completely ill. 

In favor of the finite width description one should realize that any finite width 
calculation up to a certain order includes the perturbation theory or QPA results 
up to the very same order n^. In other words, if one would expand all finite width 
terms within the order into powers of F one recovers all the QPA or perturbation 
theory terms up to that order! 

Only in the above sense current conservation, gauge invariance and other identities 
related to them (e.g. Ward identities) can be understood and expected to hold for 
a finite set of terms. These considerations are general. In particular they also apply 
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to any finite set of self-energy diagrams in the quantum case, which use full Green's 
functions with damping as discussed in the next section. 

This completes the formal discussion of the classical radiation rate and the corre- 
sponding classical expressions for the self energy of the photon. 



4 Non Equilibrium Green's Function Description 

In this section the production rate or photon self energy (2) is discussed in the 
context of non-equilibrium quantum field theory, in terms of Green's functions and 
the corresponding non-equilibrium diagrams. Thereby one has to go beyond pertur- 
bation theory which is not applicable for strongly interacting systems. The idea to 
formulate theory in terms of appropriately defined physical terms was very fruitful. 
It resulted in the development of the QPA method, sec [19,20], where the changes 
in the real part of the fermion self energy are taken into account substituting the 
free particle energy e?. by the corresponding QPA energy. This method proved to 
be very successful for the case of equilibrium matter at rather small temperature 
(T <S €f, where is the Fermi energy), since higher order corrections lead to 
additional {T/epf factors [39]. 

Since in the QPA the imaginary part of the self energy (ImE^) is supposed to be 
neghgible in the corresponding Green's functions it still suffers the same infra-red 
problems as the perturbation expansion. Any finite set of diagrams leads to infra- 
red divergences in the soft limit {q ^ 0). Even certain resummation methods as the 
hard thermal loop expansion [40,41] in QCD do not cure the problem. On the other 
hand the classical considerations of the previous section clarify, that all infra-red 
divergences disappear, if one properly accounts for the finite collision rate F. Thus, 
one has to avoid the zero-width perturbation theory or QPA and seriously account 
for the finite damping width P = — 2ImE^ of the source particles. 

Therefore the simplest and physically most meaningful step is to go to the full 
Green's functions G (full lines in diagrams) associated with the in-medium propa- 
gation for the constituents of the source. Thus one has to solve Dyson's equation 
in order to include also the damping width of the particles. This is done in sect. 
4.2. The derivation of transport schemes is summarized. Further on some physically 
meaningful resummations of the diagrams are proposed. Various graphical contri- 
butions to the self energy diagram (5) are discussed in the QPA, and in the QC 
limit. 

4-1 Vertices and Green's functions 

For a theory of fermions interacting with bosons the basic diagrammatic elements 
are the corresponding Green's functions iG (lines) and interaction vertices —iV. 



18 



For non equilibrium description it makes no sense to discuss amplitTides. Rather 
one evaluates the time-dependence of the ensemble averaged expectation values of 
observables. Such expectation values involves standard time-ordered products (T) 
for the evolution of the "ket" | ) and also anti-time ordered products (A) for the 
evolution of the "bra" ( | in the matrix elements. Both can be summarized to a 
contour-ordered product [22]. Therefore one distinguishes two types of vertices: the 
(— ) vertices with value —iV pertain to the time ordered part, and the adjoint (+) 
vertices with value -t-iV^ for the anti-time ordered section, in the here used conven- 
tion[36]. Thereby V is the real interaction vertex of the Lagrangian. Correspondingly 
two point functions, like the Green's functions, have four components, which can be 
arranged in matrix form (we reserve bold face notation for the two by two { — h} 
matrices) 



iG 



12 



( ^G^2 









T*(l)*t(2) 



T(*t(2)^(l))^ 
A^(l)^t(2)\ ^ 



(34) 



Here 1 and 2 denote the two space-time points and (...) the ensemble average. 
Upper and lower signs refer to fermion and boson Green's functions, respectively. 
The four Green's function components are obviously not independent. Rather they 
relate to the retarded and advanced ones by 



:=-i 



^(1), *^(2)]^^ 0(^1 - ^2) = G— - G-+ = G+- - G++ 



G^ +i / [^(1), ^^(2)]^\ e{t2 - h) = G— - G+- - G- 



G 



++ 



(35) 



where is the step function and [...,.. .]± denotes the fermion anti- or boson 
commutator. 

The unperturbed Green's function is resolvent of the corresponding free single 
particle Schrodinger - (non rel. fermions) or Klein-Gordon equation (rel. bosons) 



5iG?2 = (82)* GI2 = 5{l - 2)<T, , S^< 



(non rel. F) 



2m V / (3g^ 

+ A - (rel. B) 



where the subscript specifies the coordinate to differentiate. Here cr^ is the third 
Pauli matrix and 5(1 — 2) is the four-space delta function. 

With this extension to the two types of vertices — and -|- and the corresponding four 
Green's function components all standard diagrammatic Feynman rules defined for 
amplitudes can directly be generalized. These rules are then again defined with 
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respect to the real time (as in zero temperature theory)Q. For conventions and a 
detailed explanation of the diagrammatic rules used here we refer to the text book of 
Lifshitz and Pitaevskii [36] . For a given n-point function all n external vertices have 
a specified sign assignment, while one has to sum over all possible sign combinations 
at the internal vertices. 



4-2 Resummations: Dyson - and Kadanoff-Baym equations 

Since diagrammatic elements, e.g. Green's functions G and self energies S are con- 
nected at a given vertex, the same vertex type — or + appears in both functions 
and Dyson's equation can simply be written in matrix form 

Gi2 = G?2 + J d3d4 G?3S34G42 or simply as (37) 

G = G° + G°©E0G = G° + G0S©G° 

involving usual matrix algebra, which automatically provides the sum over all in- 
ternal vertex sign combinations. Here the two by two matrix S denotes the proper 
self energy of the source particlesQ. The abbreviates the space-time folding. In 
diagrams Dyson's equation becomes 




where we explicitly distinguish between fermions (label F; full straight lines for 
Gp) and bosons (label B; full wavy lines for Gb)- The four components of — iS 
are defined as the sum of all standard proper self energy diagrams like in normal 
perturbation theory, now however with definite + or — assignments at the external 
vertices, and summed over the — h signs at all internal vertices as explained above. 

Using the resolvent properties (36) of G° one can transform Dyson's equation into 
a set of integro-differential equation which in short matrix notation reads 



^ Please note that in refs. [22,23,25] the diagrammatic rules are defined with respect to 
the closed time contour, which leads to the same definition for Green's functions, however 
the off-diagonal components of the self energies used there = — S ^ and = — 
have opposite values. 

^ We reserve S for the self energy of the source particles, while 11 denotes the self energy 
of the external photon. 
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,SiGi2 = S{1 - 2)(T, + (T,j d3 S13G32 

{S2T G12 = 5(1 - 2)0-, + 1 d3 Gi3S32a, 



(39) 



The four equations involving the time changes of G~'^ and are known as the 
Kadanoff-Baym equations [23], originally derived by the imaginary time method. 
Here they are a direct consequence of the Dyson equation in matrix form. Like the 
four components of G also those of S are not independent. Their dependence can 
be determined observing that the Dyson equations for retarded or advanced Green's 
functions have to involve only retarded or advanced entities, respectively. Thus they 
completely decouple 

qR ^ Qm ^ qOR qy^Rq qR^ qA ^ ^oa ^ ^oa qj^a^ qA^ ^^q^ 

where and are the corresponding advanced and retarded self energies. Prom 
(35), (37) and (40) one therefore follows that 

= E— + E-+ , = E— + E+- , E++ + E" = - (e+- + E-+) .(41) 

The full Green's functions account for the finite damping width 

r = -2ImE^ = i (e-+ - E+-) (42) 

of the particles, which destroys the sharp relation between energy and momentum. 
Thus the spectral function 

A{x-p) := -2ImG'^(x;p) = i {G+-{x-p) - G-+{x;p)) (43) 

is no longer an on-shell 5-function by rather has a width r{x;p). This width does not 
only arises from decays (resonances) but also from collisions of the particles in dense 
matter. It is important to realize, that the resummation to full Green's functions 
(bold straight or wave lines in diagrams) reduces the set of all diagrams to a subset 
of diagrams with skeleton topology, where all self energy insertions are excluded in 
the diagrams. 

4-3 Transport equations 

In the presented picture the Wigner transforms of the off-diagonal Green's functions 
=piG'^+ and iG~^~ are Wigner densities in four space and four momentum for the 
occupied and available 'single particle states', which now have a finite width and 
therefore can be off mass shell. We will see that the corresponding non-diagonal 
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components of the proper self energy — iS are the gain and loss coefficients for the 
transport description of these Wigner densities. 

The entry to such transport equations is given by the Kadanoff-Baym equations, 
which are contained in the set (39). Subtracting the two equations for G"'^ yields 

(ia+9,J + ^) 
{-dl + dl + Ai - A2) 

~ I {^13 ^32^ + ^13^^32^ + ^13 ^32^ + ^13^^32^} 

for non rel. fermion or rel. bosons, respectively. The next step is to take the Wigner 
transformation of this equation. This involves the Wigner transformation of convolu- 
tion integrals C(l; 2) = / d3yl(l; 3)-B(3; 2) which formally can be obtained through 

C{x; k) = exp[^(a,^af - d^dl)]A{x- k)B{x; k) 

^A{x-k)B{x-k)+'-{A,B}, (45) 

where A[x'^k)^ etc. are the Wigner transforms of A, B and C and the differential 
operators act on A and B separately. The approximate expression in terms of 4- 
dimensional Poisson bracket {A, 5} = dj-Ad^B — dxAd^B defines the first order 
gradient expansion. To this order separating real and imaginary parts one obtains 
for the Wigner densities ^iG"'^ [27] 

{S{k) - E{x, k),iG-+} - {iG{x, k), E-+(x, k)} 

= k)G+-{x, k) - S+-(a;, k)G-+{x, k) (46) 

with S = (S^ + S^) /2, G = (g^ + G^) /2. 

Here S{k) — e — P/(2m) or S{k) — k"^ — is the Wigner form of the Schrodiger 
or Klein-Gordon operator (36). The first Poisson bracket on the left side gives the 
usual Vlasov part. On the right side S ^G^ and define gain and loss 

terms. This transport equation is quite general and as the original Dyson equation 
still accounts for off-mass-shell processes, part of which are contained in the peculiar 
second Poisson bracket according to Botermans and Malfliet [27]. We will see that 
the account of the finite damping width and the inclusion of higher order diagrams 
for the self energies (sect. 4.5) are essential to extend the scheme beyond classical 
transport concepts. 

Equation (46) by itself is not yet complete, since it requires a relation a relation 
between G~^ and G^ , which normally is provided by a certain physical "ansatz". 
Under near on-shell conditions, where the Wigner densities are well peaked close 
to a single on-shell energy e ~ e^-, c.f. next subsection, one can use the Kadanoff - 



^12^ 



(44) 
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Baym ansatz or apply the QPA (see Appendix A). Then the integral over e reduces 
(46) to a familiar form [25,29] 



+ vdg - {dgtj;) (9p-| n(f , p, t) 

= (±iE-+(lTn(f,p,t)) +iE+-n(f,p,t)) Zp- (47) 
= (±iE"+ - rn(f,p,i)) % 

for the on-shell particle densities of non-rel. fermions or rel. bosons 

/da f da 

— , nB{x,p,t) = i J —2eGB'^. (48) 

Here v = dj^{p) is the group velocity, and Zj^ = 1/(1— (9ReS^/(9e) for non-relativistic 
fermions and = l/(2e — dKeH^/de) for relativistic bosons is a normalization 
factor. In all expressions e is determined by the dispersion relation e = e{p), c.f. (50) 
below. 

The form of the collisional integral in (47) is convenient in many cases, e.g. to extract 
stationary solution (cf. next subsection) or for the relaxation time approximation, 
where l/FZp- is the corresponding collision time, see (A. 4). In the diffusion approxi- 
mation of eq.(47) one easily recovers eq.(14). The usual Boltzmann-hke form of the 
collisional integral is obtained from the lowest order self energy diagram which in 
QPA contributes to ImS^ (usually second order in the two body interaction), c.f. 
[25,29] and (66) in sect. 4.8. More general schemes beyond QPA arc suggested in 
sect. 4.7. The photon production rate, eq. (3), is a simple case of (47). 

4-4 Stationary and Equilibrium Properties 

For a homogeneous stationary system in Wigner {x;p) representation all Green's 
functions become space-time independent and the operation reduces to a simple 
product for the remaining 4-momentum part in the Dyson equation (37). With p = 
(e,p) we therefore drop the x-argument. Eqs. (40) can then be solved algebraically 

while G^{p) — {G^{p))* and T,^{p) — (E^(p))*. The dispersion relations are 

e + /If — ^f{^ + I^F,p), — /2m (non-rel. fermions) 

(50^ 

(e + /iB)2 = (e°^)2 + Eg(e + /iB,p), {e^^f = f + m'^ (rel bosons) 
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with corresponding free on-shell energies e^y and chemical potentials /jb and ^f- 
With r(p) = -2ImE^(p) = i - the spectral function (43) are 



^f{p) = 2 (non-rel. fermions) , 



(e + //^ - ea - ReE«(p)) + (r(p)/2)2 

Ab{p) = 2 (rel. bosons). (51) 

((e + pisr - {e^r - ReE«(p)) + (r(p)/2)^ 



They satisfy the sum rules 

de ^ , /■ , . X de 



j AB{p)2e^^l and /^^(p) 2^ = 1. (52) 



The generalization to relativistic fermions with gamma-matrices up to tensor inter- 
actions can be found in ref. [42]. For illustration and later use we give the equilib- 
rium results explicitly, which follow from the stationary condition E~+(p)G'^~(p) = 
E+~(p)G'~+(p), c.f. eq. (46), and the Kubo-Martin-Schwinger condition [15] 

E-+(p) = E+-(p)e-^/^^. (53) 



Then all the Green's functions can be expressed through either retarded or advanced 
Green's functions. From (34), (35) and (41) one then finds 

G-{p) = (1 T n,)G^(p) ± n,G^(p) , G'^ip) = ±m,A{p), (54) 
G+-{p) = -i(l T n:)A{p) , G++(p) = -(1 T n:)G\p) t n^G^'ip) , 

and the relations for the four components of the self energies 



E~ = E^ ± in,r(p) , E-+ = Tm,T{p), 

E+- =i(lTn,)r(p) , E++ = - (e--)* , (55) 

where the thermal occupations at temperature T (Fermi-Dirac or Bose-Einstein 
distributions) are 

rie = {exp[e/r] ± 1}"^ . (56) 



In the general non-equilibrium case there are no such simple relations between 
Green's functions and self energies as at equilibrium. In order to proceed one may 
simplify the problem applying so called Kadanoff-Baym ansatz (given in Appendix 
A) or using the QPA. 
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4-5 Diagrammatic Decomposition into Physical Sub-Processes 

There have been many attempts in the hterature to ehminate the redundancy in the 
definition of the four Green's function components. We do not hke to follow such 
schemes and rather prefer to keep all four components as they are, since they display 
a symmetry between the time-ordered and the anti-time ordered parts, i.e. between 
the "bra" and "ket" parts of the correlatorj^. 

In this formulation physical observables like densities, production rates, etc., are 
always given by diagrams, where the external vertices appear in conjugate pairs, i.e. 
with fixed opposite signs, just defining correlation functions. If Fourier transformed 
over space-time differences they have the properties of Wigner functions. Special 
examples are the off-diagonal components of the Green's functions and self energies, 
which define Wigner densities and gain or loss rates, respectively. Working with full 
Green's functions such correlation diagrams are given by a sum of all topologically 
distinguished skeleton diagrams with the fixed external vertices of opposite signs. 
Each diagram of given topology consists of 2" terms, where i> is the number of 
internal vertices, due to the — h summations. These extra summations make this 
approach rather non-transparent. In this section, however, we like to suggest a very 
simple classification of correlation diagrams and a reformulation of the corresponding 
sum, which amends a simple physical interpretation. 

We start with the observation, that at least in one way any correlation diagram of 
given topology and given sign assignments at all vertices can be decomposed into 
two pieces, such that each of the two sub pieces is a connected diagram which carries 
only one type of sign on all its external vertices 




The reason for such a decomposition is that then each diagram is given by a "prod- 
uct" of two sub diagrams a~ and f3~^ with the same external lines, which one may 
call amplitude - and adjoint amplitude diagrams. Here a and P denote amplitude 
diagrams including sign assignments. The adjoint of any amplitude diagram a~ 
is given by inverting the senses of all propagator lines and inverting the signs of 
all vertices; the respective values are conjugate complex to each other = a^. 

For amplitude diagrams only the external vertices have definite signs, while internal 

^ In thermal field theory the "+" vertices are often considered as ghosts. We do not like 
to support this viewpoint as these conjugate vertices are as physical as the "— " ones for 
the full correlation matrix element! 

For the construction: just deform the diagram such that all + and — vertices are placed 
left and respectively right from a vertical division line and cut along this line. Pieces which 
then become disconnected are to be reconnected to the other side until two connected sub- 
pieces remain. In case that disconnected pieces appear on both sides the result may depend 
on the order of reconnection and consequently different decompositions are possible. 
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vertices have still no sign restrictions. In analytical terms the decomposition (57) 
can be written as 



= (=Fi)G'"+(pi) . .AG^~{pm) {f3'{pi, . . . ,Pm)) a~{pi, ■ ■ ■ ,Pm) 

for space-time homogeneous cases (in general corresponding space-time foldings ap- 
pear). Here =FiG ^{pk) or iG~^ [pk) (depending on the line-sense) are the Wigner 
densities of occupied (or available states) for each external line, with off-shell 4- 
momenta pi to Pm, connecting both pieces. When a correlation diagram can be 
decomposed in more than one way according to rule (57), the picture is not unique 
and one may assign partial weights which sum to unity to the different decomposi- 
tions!^ Applying this decomposition scheme to all skeleton diagrams of the proper 
self energy leads to a decomposition in terms of physical processes with a varying 
number m of external "states" pi, . . . ,Pm besides the external photon in this case 

CO „ 

-in-+(g) = / i?"+(pi,...,pJiG-+(pi)dV---iG'+-(pJdV- (59) 

m=2 




Here the R~~^{pi, . . . ,pm) define the partial rates for each physical sub process with 
certain in- and out-states pi to pm, which all can be off-shell. In the sense of crossing 
symmetry the notion of in and out may be fixed by the line sense: incoming or 
outgoing lines. Naturally for fermions there are as many in- as out-states. It is 
important to realize that the partial rates 

i?~+(pi, . . . ,p„) = (/3~(pi, . . . ,p„)) a~{pi,...,p^) (60) 

a,l3 



arise from a restricted sum (indicated by the prime at the sum-symbol) over am- 
plitude products a times (3* such that each term arises from a given correlation 
diagram in the original sum. Note in particular that the sum (60) no /ong'er includes 
all possible interference terms of any two amplitudes! The unrestricted sum over 
all pair-products of amplitude a~ times (/9~)* is false and leads to serious incon- 
sistencies, since there are interference terms, which correspond to closed diagrams 
of non-skeleton type, which have to be omitted! The simplest case to see this is the 
normal diagram for bremsstrahlung (left) 



I 



(61) 



A simple example -for such arrcimbiguity is the following diagram which permits two 
decompositions indicated by the thin line: 
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Its absolute square leads to a diagram with self-energy insertion (61;right), which 
is not a proper skeleton diagram; rather the corresponding partial rate is already 
included in the 1-loop diagram with full Green's functions! 

Therefore, only when each term that enters a physical rate (60) originates from a 
valid decomposition of proper skeleton self energy diagrams, then this formulation 
is void of double counting. It amends a straight forward interpretation in terms of 
physical scattering processes between certain in - and out-states. These processes 
occur with partial rates R, which can be positive or negative. The warning to be 
formulated at this instance is, that if one works with full Green's functions which also 
include the finite damping width F, particular care has to be taken. The unconsidered 
account of certain Fcynman amplitudes together with finite width spectral functions 
for the particles, as sometimes done heuristically, may lead to serious inconsistencies. 

On the other hand the important point to realize at this level is that the so de- 
fined decomposition gives rise to a generalized formulation of transport theories, 
where multi-particle processes can be considered in a well defined scheme even if 
one permits for off-shell propagation. Note in particular: for a theory of fermions in- 
teracting with bosons the contribution with the fewest number of external particles 
is just three (rather than four as in the Boltzmann equation). It results from the de- 
composition a one-loop diagram and allows for one particle in and two out and vice 
versa, . Thus, in dense matter an off-shell fermion can just decay into a fermion plus 
boson or the opposite can happen, e.g. see [29,43]. For these processes it is impor- 
tant that all particles have a finite damping width in dense matter, so that creation 
and decay modes, which are forbidden from energy-momentum conservation in free 
space, may occur without principle restrictions in dense matter. 

Please note that our decomposition rules are different from the standard cutting 
rules which only apply either to the set of perturbation theory diagrams [44] or to 
the set of the quasi-particle diagrams [39]. There diagrams are cut across all H — 
lines, each cut providing an on-shell delta function from the zero width spectral 
functions. Diagrams that can be cut into more than two pieces contain a product 
of more than one delta-function on total energy conservation and therefore show a 
singular behavior which is not present in the final correlation function. Therefore 
all such "multi"-cut diagrams have to cancel out. Such arguments can no longer 
be given in the case that all spectral functions have finite widths, where the set of 
diagrams is reduced to the set of skeleton diagrams. Thus, such diagrams have also 
to be considered in a description with full Green's functions! 

4-6 Three and Four Point Functions 

In principle one can stick to the above picture since all infra-red divergences disap- 
pear for all diagrams due to the finite propagation times of all Green's functions. 
However one may have to consider still quite a numerous amount of diagrams in or- 
der to achieve meaningful results. For instance on the QP level we expect, that the 
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proper in-medium current appears at both external vertices in a symmetric fashion. 
Thus for the convective currents one expects 

f ^e.' = e±-.m^e(±-4^ ^R.E«) / (l - |ReE«) , (62) 

where e(p) — e^-\- ReE^(e(p),p) defines the QP energy momentum relation, see 
eq. (A. 3) in Appendix A. In diagrammatic terms this can be achieved by certain 
partial resummation that lead to vertex corrections. Thereby we shall not consider 
an immediate resummation to the exact full vertex, since this would amount to solve 
the whole problem. Rather we like to stay to a picture where in certain limits like 
the QPA and QC limit, piece by piece an interpretation in physically meaningful 
terms can be given. 

The considerations in the preceding subsection assigned a particular role to the full 
— h and H — Green's functions as Wigner densities. This suggests to apply further 
resummations and to extend the ideas put forward in ref. [39] in the context of quasi- 
particles now to particles with finite width. Namely, one likes to gather diagram 
pieces that are void of the Wigner densities G^^ or , both for fcrmions and 
bosons. That is, one likes to rcsum sub-pieces of skeleton diagrams with given number 
and type of external vertices (3 point or 4 point functions, for example) where all 
internal and external vertices have only one definite sign value. The {— } diagrams 
then contain only time-ordered full Green's functions G and therefore represent 
a straight forward generalization of the standard zero temperature Feynman 3 or 
4 point functions now including the full self energies. The {-|-} diagrams are just 
the adjoined expressions. This way one can define 4 point functions (in-medium 
interactions) 




Here the thick wavy lines relate to the corresponding G^^-exact boson propaga- 
tors or two-body potentials in non-rclativistic theories with potentials. Since only 
like sign vertices are permitted, no G^^ and G^~ lines appear in these functions. 
Such resummed expressions have been proven useful in the low temperature QPA 
to define in-medium interactions and effective vertices and they appear also quite 
meaningful in the limit of low densities, as in the classical limit for example. Vari- 
ous approximation levels are possible for the 4-point functions; a detailed discussion 
would be beyond the scope of this presentation. We mention just a few possibilities: 

a) a ladder summation in the s-channel (horizontal in diagram (63)) for the particle- 
particle (p-p) and particle-hole {p-h) channels generalizes the Bruckner G-matrix 
to non-equilibrium; 

b) in many practical cases (e.g. Landau - Migdal's Fermi-liquid theory [19,20]) the 
4-point functions are approximated by 2-point approximants in the t-channel (ver- 
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tical in (63)). Then RPA-type resummations are possible [39,10], which iterate 

" " and "+ +" loops in the t-channel; details are given in Appendix B; where 

also the scheme of the bosonization of the interaction is presented, cf. [10,48]. 
d) the ultimate could be a crossing and exchange symmetric form; in this case one 
relies on suitable parameterizations. 

These like-sign effective interactions generalize the two-body scattering matrix in 
matter to non-equilibrium. Thereby one does not only account for the change of the 
fermionic occupations (as already considered in the literature) but also includes the 
damping of the fermions. In this respect it would be interesting to see, how bound 
states (e.g. the dcutcron [49] or the J/^) change their properties in dense matter 
(Mott transition) also due to the damping widths. For consistency these like-sign 
effective interactions then also enter the definition of the in-medium vertices (3 point 
functions) defined as 



(64) 



4-7 Key diagrams 



For simphcity we confine the discussion to the case where the 'external' photon just 
couples to fermions. Any generalization to other types of particles (external and 
source internal) is straight forward. Due to the above considerations from now on the 
remaining diagrams include the following elements: full fermion Green's functions, 
like-sign 4-point interactions and the corresponding vertices. Please notice that this 
reduces the set of diagrams even further! In particular not all sign combinations are 
permitted any longer since some of them are already included in the resummed 3 or 
4 point functions. 

All photon self-energy diagrams can be build up by iterative four point insertions. 
Thus, the set of diagrams for reduces to 

This set of "key"- diagrams is important for all subsequent considerations and there- 
fore deserves further comments. 

(1) Each diagram in (65) represents already a whole class of pcrturbativc diagrams 
of any order in the interaction strength and in the number of loops. The most 
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essential term is the one-loop diagramP^, which is positive definite, and cor- 
responds to the first term of the classical Langevin result for Hd in (24) as 
we shall show later. The other diagrams represent interference terms due to 
rescattering. 

(2) Compared to conventional diagrams, vertex corrections can appear on both 
sides of one loop as they are separated by {H — } lines (see example given in 
Appendix B). 

(3) Note that the restriction to like sign vertices for the resummations (64) and 
(63) are defined with respect to skeleton diagrams in terms of full Green's 
functions. "Opened" to perturbative diagrams with thin lines, these can 
still contain alternative signs, since Dyson's equation (37) includes all signs in 
the intermediate summations! 

(4) In some simplified representations (being often used) the 4-point functions be- 
have like intermediate bosons (e.g. phonons), c.f. Appendix B. 

(5) For particle propagation in an external field, e.g. infinitely heavy scattering 
centers, only the one-loop diagram remains, since the one deals with a gen- 
uine one-body problem. However, extra complications arise, since translation 
invariance is generally broken and the Green's functions then also depend on 

X. 



4-8 Decomposition of Closed Diagrams into Feynman Amplitudes in the QPA 

The QPA is quite commonly used concept originally derived for Fermi liquids at 
low temperatures (Landau-Migdal, see [19,20]). There one deals with on-mass-shell 
fermions in matter (quasi-particles) described by the pole part of the Green's func- 
tions, i.e. one assumes that ImS^ — > in the Green's function Gp. Then with the 
help of some phenomenologically introduced interaction (particle-hole irreducible) 
one calculates the values ReS^ and ImS^ which now depend on quasi-particle prop- 
erties. Since in QPA the finite width contributions have to appear in higher order 
through corresponding ImS-insertions the whole set of QPA diagrams defining the 
full — in""*" is by far larger, than set (65). 

The QPA has considerable computational advantages as Wigner densities (" — h" 
and "H — " lines) become energy 5-functions, and the particle occupations can be 
considered to depend on momentum only rather than on the energy variablep*^. 
Formally the energy integrals in eq.(5), (65) can be eliminated, in diagrammatic 
terms just cutting the corresponding "- +" and "H — " lines [39]. This way one 
establishes a correspondence between correlation diagrams (65) and usual Feynman 
amplitudes in terms of QPA asymptotic states and QPA Green's functions. Thus 

In perturbation theory or QPA the corresponding one-loop diagram usually vanishes 
for on-shell photons due to conservation laws. Here however, with full Green's and vertex 
functions it represents a series of perturbative diagrams as the reader can easily imagine. 

The later approximation is also often used beyond the scope of the QPA and is then 
known as Kadanoff-Baym ansatz, see [23] and Appendix A, c.f. [12] 
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the QPA allows a transparent interpretation of correlation diagrams. 



For the dynamics of the fermion transport the first QPA diagram that contributes 
to the gain (loss) term 



/^E-+(Pi)G--(Pi)-{,Q_}^p^(l-nO (66) 
j d'^p2 . . . d V 



QPA 

|2 

TT 5'^(pi +P2-P3 -P4)(l - ni){l - n2)mn4 



leads to the standard Boltzmann collision term with corresponding occupation and 
Pauli suppression factors for the in and out states. Here and below the full blocks 
denote the effective two-fermion interactions, and thick fermion lines denote the 
QPA states or Green's functions. 

For the here studied photon rates we discuss in detail the correlations diagrams on 
the right side of (65) with consecutive numbers 1 to 6. Thereby diagrams 1, 2, 4 
and 5 describe the bremsstrahlung related to a single in-medium scattering of two 
fermionic quasi-particles and can be symbolically expressed as Feynman amplitude 
(67a) 

(a) '* : (b) -- V\> -"-- (67) 



(c) J :■ '* : id) 



The full circle denotes the effective vertex. One should bear in mind that the photon 
may couple to any of the external fermion legs and all exchange combinations are 
possible. The one-loop diagram in (65) is particular, since its QP approximant van- 
ishes for real or time-like photons. However the full one-loop includes QPA graphs 
of the type (67b), which survive to the same order in F/e as the other diagrams [39]. 
In fact it is positive definite and corresponds to the absolute square of the amplitude 
(67a)), c.f. (61). The other diagrams 2, 4 and 5 of (65) describe the interference of 
amplitude (67a) either with those where the photon couples to another leg or with 
one of the exchange diagrams. Thereby for neutral interactions diagram (65:2) is 
more important than diagram 4 , while this behavior reverses for charge exchange 
interactions (the latter is also important for gluon radiation from quarks in QCD 
transport due to color exchange interactions). Diagrams like 3 describe the interfer- 
ence terms due to further rescatterings of the source fermion with others. According 
to our rules the diagram (65:3) corresponds to a two-body collision process and de- 
scribes the interference of amplitude (67c) with that one where the photon couples 
to the initial leg. Diagram (65:6) describes the photon production from intermediate 
states and is given by Feynman graph (67d). In the soft photons limit {ujg ^ e^) this 
diagram (67d) gives a smaller contribution to the photon production rate than the 
diagram (67a) in QPA, where the normal bremsstrahlung contribution diverges like 
l/uiq compared to the l/e^-value typical for the coupling to intermediate fermion 
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lines [45]. However in some specific cases the process (67d) might be very important 
even in the soft hmit. This is indeed the case for so called modified URCA process 
nn — >■ npeV which is of prime importance in the problem of neutrino radiation from 
the dense neutron star interior, see [43]. 

Some of the diagrams, (c.f. the graph shown in footnote 7) which are not presented 
explicitly in eq.(65) give more than two pieces, if being cut, so they do not reduce 
to the Feynman amplitudes. 

For the validity of the QPA one normally assumes that F <^ e, where e is an average 
particle kinetic energy (~ T for equilibrium matter). With F ~ Tr^T^/ej? for Fermi 
liquids, c.f. (A. 5), the QPA constitutes a consistent scheme for all thermal excitations 
Ae T <^ ep. However with the application of transport models to higher energies 
this concept has been taken over to a regime where its validity can no longer easily 
be justified. Moreover, our considerations show that the condition F ^ e is not 
at all sufficient. Rather one has to demand that also ^ F in the QPA, since 
finally energy differences of order u appearP^. In particular, the remaining series 
of QPA-diagrams is no longer convergent unless u > T, since arbitrary powers in 
r/uj appear, and there is no hope to ever recover a reliable result by a finite number 
of QPA-diagrams for the production of soft quanta! With full Green's functions, 
however, one obtains a description that uniformly covers both the soft {uj ^ F) and 
the hard {uj ^ F) regime. 



5 The Quasi-classical Limit (QC) 

In this section we like to discuss, which class of diagrams remains in the quasi- 
classical limit and how this is to be interpreted. 

The QC limit requires that 

i) all occupations of the source particles are small ^ 1) implying a 

Boltzmann gas with fi ^ —T and that 
ii) all inverse length or time-scales times h are small compared to the typical 
momentum and energy scales of the source systems. 

In particular this implies hu, <C e, and a collision rate F = h/Tcow -C e, where 
e is a typical particle kinetic energy (~ T for equilibrium matter). To be precise: uj 
and q of the produced particle are sensitive to the space-time structure of the source, 
while they are negligible as far as energy and momentum balances are concerned. The 
latter fact permits to prove the Kadanoff-Baym ansatz in this case (see Appendix A 
and discussion of eq. (84) below) which considers the occupations of the source 
particles to dependent only on momentum n^- = n^^^^^ but no longer on energy e. 

This statement is particular, since one compares the photon energy to with the damp- 
ing width of the source particles F, while the damping rate of the photon itself 7 ~ 
|n ^1 /{2ujn^) < A-ire"^ po/{Tm) ~ 20MeV^/r for nuclear matter can be quite small! 
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Also wc assume that F will not depend on time in between subsequent collisions 

(|r| ^ 1/r). 

We note in particular that for bosons with chemical potential hb — 0, hke the 
produced photon, the equilibrium occupations will be large, ub ~ T/cu ^ 1! This 
fact is of no further relevance, if one excludes internal photon lines in the proper 
correlation functions (2), n~+ and n+~ as we do. 

5.1 Time Structure of Green's Functions and Loops in QC Limit 

For fermion Green's functions one has the following simplifications 

Gj^{p)^mpA{p), G'+-(p)--i(l-np)^(p), (68) 

while at large temperatures T the particle occupations are given by 

np ~ exp[-(ep- - /x^)/T] < 1. (69) 

The correspondence between the diagrammatic expansion (65) and classical limit 
of sect. 3 becomes more transparent if one uses the mixed r — p representation for 
the Green's functions, where r is the time difference between the two space-time 
points. Then from the definition of and Green's functions (see eq. (34)) 
one immediately finds for fermions 

G'^+(t,p) ~ inp-exp[-|FT|/2 - ie^], (70) 

G+-{t,p) ~ -i(l - np) exp[-|FT|/2 - ie^], (71) 

while Gp~ = (1 — np)Gp — n^Gp and Gp'^ = — (1 — np)Gp + npGp are essentially 
retarded and advanced, respectively. 

A further simplification comes from the time behaviour of fermion-loops, c.f (B.7), 
which mediate classical energy and momentum transfers. The corresponding time 
scales 1/e or l/{vp) are very short on the damping scale 1/F, so that such loop 
insertions become instantaneous. This is the reason, why one recovers a Markovian 
description for the motion of the source in the classical limit. With these simplifica- 
tions we now calculate the diagrams (65). 

5.2 Self Energy Diagrams in the QC Approximation 

In the mixed r — q representation the one loop diagram is given by 
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-ino-+(r,g-)= I ^^L^yMy^G'-+(r,pOG+-(r,pl)(27r)-^<5(pl-g-pl) 

^ / W)'^^'^'^''^^^ ^/^^^^ ~ ?72))e-l^^le-'(«~)^ (72) 



if \qv\ -C T. Here ~ yA^j^{p) defines the in- medium photon - fermion vertex 
in the classical limit following eq. (62), while p ^ {pi +p2)/2. Apart from the q- 
dependent oscillations the time structure of this diagram is given by an exponential 
decay: e"''""^' which leads to 



j7r^n{p + q/2){l - n{p- q/2))- 



(27r)3'"^" ' ^^(a;_gy))2 + r2 

for the spatial components of n~+. This expression is identical to the n — term of 
the classical Langevin result (25). 

The classical Langevin example (sect. 3) considers the propagation of a single charge 
(say a proton) in neutral matter (e.g neutrons). Therefore for this case only diagrams 
occur, where both photon vertices attach to the same proton line. Also, of course, 
no direct proton-proton interactions occur. In the following we like to show that 
diagrams of the type 



- +(74.) 



with n { — h} scattering interactions (intermediate particle-hole { — h} neutron loops) 
correspond to the n-th term in the Langevin result (24). To demonstrate this we look 
into the time structure of such a diagram, and assign times and r to the external 
— and -|- vertices, while the — and + interactions arc taken at tj~ to t~ and to t^, 
respectively. In the classical limit G is retarded, while G'^'^ is advanced (see eqs. 
(35), (70) and (71)), such that both time sequences have the same time ordering: 
< ti < . . . < < T and < tf < . . . < t:^ < T (all inequahties reverse, if all 
line senses are reversed) . Thus the r-dependence of the modulus of these diagrams 
gives e~l'"^'. A second simplification emerges from the fact that the H — loop inter- 
action insertions mediate classical momentum transfers \pn — Pn+i \ which are large 
compared to ^F. Therefore the time structure of these loops becomes very short 
on the scale l/F and therefore merge (5-functions: S{t~ — t^). With tn = t~ = t+ 
diagram (74) then no longer depends on the intermediate times tn apart from the 
ordering condition, and therefore results in a factor |rr|"/n! With hq — also the 
corresponding momenta are pair-wise identical, and the remaining momentum inte- 
grations just serve to define the correlation between Vm and Vm+n after n scattering. 
Thus 

+ + + 

-- {([)([)••• 4)} -"- (75) 
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' -in-+(r, g = 0) 47reVo ^eH^-l 

^ -in-+(a;, g = 0) ~ 47reVo (^m^'m+n)^ 2r"Re y_.^)„+, 

where the resulting proportionality to the proton density po results from — h and 
H — Green's functions next to the external vertices. Here of course we have silently 
assumed a consistency condition to be fulfilled: namely that the interaction loops 
which one takes into account also consistently define the damping width F of the 
source particles! In particular the extra powers in occupations coming from the H — 
loops are contained in F, and therefore do no longer explicitly appear in the final 
result. This proves that in the classical limit these diagrams reproduce the terms of 
the classical Langevin series (24). 



5.3 Hierarchy of the QC expansion 



The Langevin diagrams (75) have the following properties: 

a) external vertices: the external photon couples to the same fermion line; 

b) topology: the diagrams are planar, i.e. no crossing of lines occur. 

c) — h sign topology: if one cuts them at all — h lines they decompose precisely 
into two pieces; 

d) value: apart from the velocity correlations they all give the same contribution 
to the soft photon point (cu, q) = 0.; 

Using the equal time properties of classical — h- loop interaction insertions one finds 
that 

+ oo + - oo 

i) = j Fe-'^^dr = 1, while (i = j nFe'^^dr = n (76) 

"*! T 



at g = 0, where n is the proton occupation. The left case iteratively enters the 
Langevin diagrams and one compiles factors of unity for each time folding, since the 
loop insertions of order F are compensated by the time integration over the Green's 
function product G This proves d). 

What remains to be shown is that for the classical problems discussed in sect. 3 all 
other key diagrams are disfavored by extra occupation factors n or F/e which both 
are small compared to unity. 

To a): Since in the classical problems of sect. 3 only a single charged particle (say 
a proton) in neutral matter is considered, all diagrams with more than one proton 
line do not occur in this case. To b): Non-planar diagrams, where interactions cross, 
violate classical time-ordering. Either the two interaction times are interlocked on 
a time scale 1/e and therefore lead to a penalty factor F/e or restoring the time- 
ordering one obtains a Z-shape fermion line, which gives an extra occupation factor 
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n. To c): Diagrams that can be cut into more than two pieces can be obtained on 
various ways: i) by any — h- loop insertion which switches signs on the outer proton 
hues, thus using the right insertion in (76); here one obtains an extra hindrance 

factor n. ii) extra insertions of or ++ blocks linking the outer proton lines; 

since there are no direct proton-proton interactions these effective interactions are 

mediated by the surrounding neutron matter thus containing neutron loops, 

which also leads to an additional factor n, besides another n factor due to the — h 
proton Green's function. 



6 General Quantum Consideration for Hot and Dense Matter 

The considerations above show the following. A proper treatment of an entirely 
classical problem, namely the coupling of a classical source to a wave (electromag- 
netic field), on the level of quantum many-body theory requires technics, that even 
nowadays are still non-standard, i.e. beyond perturbation theory or QPA. While the 
classical problem can be solved quite conveniently and simple with no problems on 
the infra red side, the corresponding quantum description requires an appropriate 
account of the finite damping width F of the source particles. The most natural ap- 
proach in our mind is the real-time Green's function technic, which however requires 
partial rcsummation, such that the finite width is included already on the one-body 
Green's function level. 

In this section we analyze the production rate from hot and dense matter in the 
quantum case in terms of non-equilibrium Green's functions. In order to provide 
some analytical results which easily can be discussed in different limiting cases, 
we employ the following approximation for the full retarded Green's function. We 
assume to be given by a simple pole approximation with constant residue 

^ e + i^F- e{p) + ir/2 

where e{p) = p^/{2m*p), m*p being an effective fermion mass. The width T is assumed 
to be independent of e and p. Explicit results will be given for the one - and three 
loop case (the first two diagrams in (65)). 

6.1 Qualitative expectations 

Compared to the classical results we expect the following changes: 

a) one looses the classical hierarchy of diagrams, such that many more diagrams 
contribute in the quantum case; for specific couplings, however, some diagrams 
are disfavored or drop due to selection or suppression rules, e.g. non-planar 
diagrams in SU(n) coupling; 
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(77) 



b) the radiated quantum carries finite momentum and energy (which vanish in the 
classical limit), such that additional recoil corrections and phase-space factors 
~ g-'^/r appear; the latter is important, since it cures the classical ultra violet 
catastrophe, where the intensity spectrum is white, leading to a divergence of 
the radiated energy in the classical case; 

c) the occupations are no longer of Boltzmann type with n <C 1, so that Pauli 
suppression and Bose-Einstein enhancement effects are significant; 

d) the duration time of binary collisions ~ l/e^ mediated by interaction loop 
insertions (B.7) in the correlation diagrams, are no longer negligible compare 
to 1/r as in the classical case, so that non-markovian memory effects become 
important [46,47]. 

Points (b) and (c) can be clarified using the exact relation between thermal fermion 
and boson occupations 

n^(e + a;/2)(l - n^(e - uj/2)) = (n^(e - uj/2) - n^(e + a;/2))n^(a;) (78) 
which simplifies in the following limits to 



{e + uj /2){l -n^{e-uj/ 2)) ^-^y{e)n^{uj)uj for uj (79) 

ftin^(e)n^(cu)cu/T for cu<T,n^<l. 

At low temperatures only states close to the Fermi-surface contribute. The last 
approximate relation suggests that relative to the classical results of sect. 3 an 
additional phase-space suppression factor n^{uj)uj/T appears in the quantum case, 
which accounts for the finite energy uo carried away by the quantum. 



6.2 Contribution of One-Loop Diagram with Full Fermion Propagators 

We first consider the one-loop diagram of (65) with the full fermion propagators 




= -in-+ = -iV^VAo^ (80) 



where ^ denotes the bare loop without vertices. For simplicity we will neglect ver- 
tex corrections. The later can be trivially included in Landau-Migdal approximation 
(e.g. see eq. (72), Appendix B and refs. [10,48]). 

With the help of relation (78) and the equilibrium form of the Green's functions 
(54) the bare loop reads 

- i^-+ = I ^ - ^ « - (81) 

J (27r)4 (e + /x^ - + {V/2f (e + /x^ - e-^,^ + + {V/2f ' ^ ^ 
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This expression can be evaluated in closed form in different limits. We first analyze 
the QPA ( i.e. F <^ u,kvF,T in case T -C or F <^ oj.kvr ~ k\jT /m*p, for 
T ^ ep)-, which can be found in the literature. In this limit one recognizes two 
energy ^-functions in eq. (81) (see approximation formula (A. 2) of Appendix A), 
which together with momentum conservation and exact relativistic kinematic can 
only be fulfilled for space-like {uj,k). Following ref. [10] one obtains for the QPA 
loop in various limits 



-i Aq^{uj, k) 



QPA 



* \2rp 



FJ - „B ^^ exp{K) + 1 

2^1^''^ exp(K) + exp(-u;/r) 



for u > \k\ 
for <^ l/cl 



(m 



* \2 



T 



~^Tfk 
(m 



-n, 



B 



for /t ^ — 1, 



2nk ^^P(~^) exp(-t^/r) for k > 1, 



id < \k\ 
uj < \k\, 



(82) 



where k = {m*pv'^ /2 — fip)/T, with v = — k'^ / {2m*p))/\k\. 

Here non-relativistic kinematicsp^ has been used, where v is the recoil corrected 
fermion velocity that essentially contributes to the loop and the condition u <^ \k\ 
assures |f | <C 1. The simplified expression for k <^ — 1, realized for T -C ei? ~ /i^? and 
UJ < kvp, is quite frequently used for space-hke interaction loops in low temperature 
Fermi systems as in Landau's Fermi liquid theory. However, it shows a singular 
behavior oc l/{kuj) in the small <^ = |fc| -C T limit, which is a generic defect 
of the QPA. The Boltzmann limit fip —T leads to the k ^ 1 case, which apart 
from recoil correction and the extra quantum phase-space factor e~^^^ coincides 
with that for classical diffusion result (20) in sect. 3. In summary, for the one loop 
term the QPA leads to meaningful results only for large space-like u, q (hard thermal 
loops). 

In the general case only the integral in eq. (81) over the angle pk can be performed 



oo oo 



n, 



de 



pdpm*p 



F 



X 



arctan 



— arctan 



(27r)3fc {e + fiF- e^f + (F/2)2 
e + /iF + - Cp - ejt + pk/m*p 

e + /iF + - Cp - efe - pk/ m*p 



n 



(83) 



(F/2) 



in closed form. To proceed further we consider the case of small spatial momenta 
k = \k\ such that pk <^ m^F. For the remaining two-dimensional integral 



As a defect of the non-relativistic approximation the result does not exactly vanish for 
LO > \k\ but rather leads to terms of the order exp{—m*p/T) or less. 
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, l_ f ^TTt jpj Fi f -I / W \ f ^ / \ 



oo 

>l/2_ 



— oo 

r 

X 1 ^ 7^77^ for pk -C m%V (84) 

(e + /XF-ep-efe + u;)2 + (r/2)2 ^ ^ ^ ^ 

one realizes that the ep-integration over the product of Lorentz functions gives 



A^T^TTjr, e + /.^»r (85) 



(a;- A;2/(2m»)2 + r2 



while it is essentially zero for e + /xf -C — F. For correspondingly small values of F 
one therefore obtains 



T ( \ 3 /2 

- - 27r2[(u;-A:2/2m^)2 + r2] J "^"f ^''^.-mf " ^e,+c.-^J, (86) 



where we have replaced the remaining e variable by e^. The very same form emerges, 
if one formally replaces 

nf - nf+^ ^ - (87) 

in (84) and first integrates over e. This approximation (87) corresponds to the 
Kadanoff-Baym ansatz (see Appendix A). 



Equation (86) vahd for pk <^ m*pV can be evaluated in two limits 

" T "(cu-A;2/(2m>))2 + F2 

where essential contributions arise from momenta ~ pi? for T <^ or around 



rrfpPFTlTi'^ for F,T<eF~A*F , , 

(88) 

Pf for T»eir,a;,F 



Pt ~ ym*pT for T > ep. Here p^? is the density of the charged fermions. Compared 
to the QPA which is zero for time-like momenta, this result is finite and of order 
l/F in the soft limit. It agrees with the classical result (25) and the corresponding 
QC limit (73) besides recoil and the quantum phase-space corrections. 

Starting from the QPA for the fermion Green's functions one usually attempts to 
restore a dependence on the non-zero fermion width for the boson self-energy 
by means of the analytical continuation a; ^ a; -|- iF, where 

.R_ f d^P - ^p+k 



-^''J {2nrep-e,^, + u; + ir- ^^^^ 

In refs. [50,10,51] such a procedure has been used in order to account for the finite 
A-isobar width in the pion self energy. We see that the value Im^^ given by eq. (89) 
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and that given by (86) coincide only if n-^-^ ~ nep+oj-^ip, i-e. applying Kadanoff- 
Baym ansatz. 

The thus discussed one-loop diagram can also be used for the intermediate t-channel 
interaction loops Q~'^ occuring in higher order diagrams. There the typical values 
of parameters are 

u! < ep , k > pf , p > Pf ■ (90) 

Therefore at least for F ^ ei? one has large space- like momenta with pk ^ Fm^. 
Then from eq. (83) one obtains 

- '^^^ = - / ^^^^ - + arctan(^)] , (91) 

— oo 

which merges the QPA expressions in (82) for large space-like momenta in the limit 
F <C ej? both at low temperatures T <^ ej? ~ /xj? (k <C —1) and in the Boltzmann 
limit («>!). 

Thus, for the intermediate t-channel interaction loops Q~~^ one can safely use the 
QPA (F — > 0), which even accounts for higher order correction to this loop, as one 
sees from the corresponding hmit of the diffusion result (20). Only for soft loops 
(e.g. as in the s-channel) the QPA is ill defined. 

Comparing the one-loop result at non-zero F (88) with the first non-zero diagram 
in the QPA (F = in the fermion Green's functions) 




at small momentum k one determines a correction factor 

W^^;;^, (93) 

which cures the defect of the QPA for soft lo. This factor comphes with the re- 
placement a; ^ a; -|- iF. A similar factor has been observed in the diffusion result, 
where however the macroscopic relaxation rate F^; enters, due to the resummation of 
all rescattering processes. Other factors between eq. (88) and in the corresponding 
QPA Feynman diagram (c.f. ref. [45]) become identical, if one explicitly calculates 
the width F to that order, see also [50] . 
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6.3 Higher Order Diagrams with Full Fermion Propagators 



Along similar routes (details are given in Appendix C) the correction factors for the 
higher order diagrams can be derived. Here we just quote the results for the next 
lowest order diagrams 



■■*<$>"-^''-'{--*<i>"- 



■'00" 



with Co (a;) from (93) and 



ui" 



QPA 



QPA 



(94) 



(95) 



(a;2 + r2)2 



(96) 



The total radiation rate is obtained from all diagrams in (65) . 

The set of QPA diagrams in (92), (94) and (95) are just those that determine the IQF 
scattering rate (28) including the exchange diagram for a source of fermions (95). 
The latter drops in the classical limit, where the two other ones yield the n = and 
n = 1 terms of the classical Langevin result (24). Thereby the damping correction 
factors Co and Ci in the quantum case are the same as classically derived. For 
r <^ a; they tend to unity and the production rate coincides with the QPA results 
as obtained in ref . [45] , while there is a substantial suppression at small frequencies 

In the corresponding QC limit all the diagrams of type (74) with an arbitrary number 
of — h A'"A'"-interaction insertions can be summed up leading to the diffusion result 
(18) in sect. 3.1. For small momenta q this leads to a suppression factor of the form 

There is hope that even in the quantum case some higher order diagrams can also 
be resummed and that qualitatively a similar suppression factor emerges like for the 
diffusion result. 



7 Discussion and Perspectives 



We investigated the production of particles from the collision dynamics of dense mat- 
ter at the example of photon production. Thereby the source of charged particles 
was described in two ways, a) as a classical system governed by classical transport 
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equations and b) as a quantum system in terms of a real-time non-equilibrium field 
theory formulation. The central quantity is the current-current correlation function 
which relates to the imaginary part of the proper self energy of the produced par- 
ticle. Under quite general assumptions this correlation function governs the local 
production and absorption rates in the matter. Since the here discussed features are 
of kinematical origin, relating space-time scales to the corresponding momentum 
and energy scales, all conclusions drawn in this study are general and therefore also 
apply to the in-medium production and absorption rates of any kind of particle. 

The problem could be quite naturally formulated and solved in the classical de- 
scriptions by means of a macroscopic transport and a microscopic Langevin process. 
These studies showed that spectrum of produced particles is essentially governed by 
one macroscopic scale, the relaxation rate r^; of the source. For frequencies u) of the 
produced quantum which are large compared to the spectrum can be described 
by the incoherent quasi-free scattering approximation (IQF) used in most of the 
transport models. Higher order corrections help to improve the result. Once, how- 
ever, u; < this quantum is "soft". It can no longer resolve the individual collisions 
in time and therefore the IQF picture fails and produces a false infra-red divergence 
in the rate. Rather, the correct rate is regular and differs from the IQF result by a 
suppression factor 



This soft part of the spectrum is genuine non-perturbative. The essential features 
are summarized in figs. 1 and 2 of sect. 3. For relativistic sources a second scale 
comes in, once the wave number \q\ of the produced particle exceeds the value of 



Vx/ y {v'^) due to the increased spatial resolution as shown by the closed form results 
of the diffusion model. 

On the quantum level new scales come in since now uj and q also correspond to 
the energy and momentum of the photon which have to be compared with the 
characteristic energy and momentum scales of the source as given by temperature 
T and chemical potential or the Fermi energy ep- Also the occupations can be- 
come degenerate and Pauli suppression or Bose enhancement effects are important. 
For the general formulation one has to leave theoretical schemes that are based on 
the concept of asymptotic states like perturbation theory or quasi-particle approx- 
imation. The proper frame is the real-time non- equilibrium field theory, where the 
proper self-energy can be formulated in terms of closed correlation diagrams with 
general propagators. Thereby the resummation of Dyson's equation to full propaga- 
tors which also include the imaginary parts of the self energy and therefore account 
for the damping of the source particles is the essential step to cure the infra-red 
problem. Only this way one comes to a convergent scheme. After this resumma- 
tion the corresponding set of diagrams is then reduced to diagrams with skeleton 
topology. Using the Kcldysh — h notations we have addressed a particular role to 
all — h and H — fines as 8-dimensional Wigner densities of occupied and available 



Co{uj) 



a;2 + r2 ■ 



(97) 
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" states" . This motivated further resummations which define in-medium interactions 
and vertex corrections. The resulting set of diagrams can then be discussed in detail 
in these physical terms. A decomposition of these correlation diagrams in terms of 
the interference of two "amplitude" diagrams is suggested, where Wigner densities 
enter as in- and out-states. This permits a transparent physical interpretation of 
the correlations diagrams, which may be used to formulate multi-particle collision 
processes in matter. They can serve as input for a generalized transport description, 
which ultimately includes the off-shell propagation of particles and therefore unifies 
resonances which have a width already in vacuum with all other particles in the 
dense matter, which acquire a damping width due to collisions. 

Once one has the full propagators and the in-medium interactions it is in principle 
straight forward calculate the diagrams. However both, the computational effort to 
calculate a single diagram and the number of diagrams, are increasing dramatically 
with the loop order, such that in practice only lowest order loop diagrams can be 
considered in the full quantum case. In certain limits some diagrams drop out. In 
particular we could show, that in the classical limit of the quantum description 
only a special set of diagrams survive, which could be associated with the multiple 
collision terms of the classical random Langevin process. Comparing the lowest order 
loop diagrams in various limits to the corresponding QPA diagrams one realizes that 
also here correction factors similar to (97) appear. Now the characteristic scale is 
the damping width F of the source particles. Accounting for higher order diagrams 
one concludes that also in the quantum case the relaxation rate r^; is the relevant 
scale which decides between soft and hard photons. Thus for applications one has 
to compare the typical energies of the produced particles with the typical relaxation 
rates of the source system. 

Our considerations are of particular importance for the theoretical description of 
nucleus-nucleus collisions at intermediate to relativistic energies. With tempera- 
tures T in the range of 30 to 100 MeV for dense nuclear matter, up to 200 MeV for 
hadronic matter and beyond 150 MeV for the quark gluon plasma or parton phase 
most of the kinetic models that are used infer collision rates T for the constituents, 
which during the high density phase can reach the system's temperature, F ^ T. 
Such estimates make the use of on-shell concepts already rather questionable. The 
particles uncertainty in energy is comparable with the mean kinetic energy! In par- 
ticular the bulk production and absorption rates of all particles with masses less than 
T, if calculated in standard IQF approximation, are seriously subjected to the here 
discussed effect. Therefore the corresponding quenching factors (97) should sensi- 
tively affect the production rates of quark pairs and gluons during the plasma phase, 
of low energy pions during during hadronization and real and virtual photons with 
correspondingly low energies. Since our discussion was restricted to the production 
in dense matter, for the particular case of photon production in nuclear collisions 
one has to consider in addition the radiation caused by the incoming charged ions 
and outgoing charged fragments. Due to Low's theorem [55] the latter give rise to an 
infra-red divergent ^ 1/ui component which interferes with the one discussed here. 
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In astrophysics neutrinos produced from neutron stars or during super nova collapse 
have an absorption mean free path which is long compared to the size of the radiating 
system (for review see [7,8,10]). Thus the production rates cannot be estimated by 
black-body radiation. Rather the microscopic rates are relevant. The mean kinetic 
energies per neutrino or uP-pedr are about ~ ST or ~ 6T, respectively [52,53]. With 
mean collision rates of the order of F ~ Ti'^T'^/ep <C T [54], c.f. eq. (A. 5), the 
production rates can safely be estimated in IQF-approximation for relatively cold 
neutron stars T < 1.5 to 2 MeV. Already around T ~ 5 MeV the quenching factor 
(97) is significant (0.3) and it may become even smaller during super nova collapse. 
Then the temperatures can raise to T ~ 10 — 30 MeV such that F ~ T and the here 
discussed suppression effects are relevant for the corresponding neutrino emissivity. 

Appendix 



A KadanofF— Baym Ansatz and QPA 

In the general non-equilibrium case one has no simple relations between Green's 
functions as in equilibrium. In order to proceed nevertheless one often uses the so 
called Kadanoff-Baym ansatz [23]. For Fermions it reads 

G+- = 2i(l - nj__^ JImGP, = -2inJ._^^lmGP, (A.l) 

where n^^_^^ are the fermion occupations which now depend on p through the on- 
shell dispersion relation (50) rather than on e. One should note that the Kadanoff- 
Baym ansatz does not directly follow from the properties of the and 
functions, rather it has been introduced in order to recover the Boltzmann limit. 
The correctness of this ansatz has only been proven in the QPA, see [12]. Eqs. (A.l) 
complies with the definition of the particle densities (48), as can be seen by the sum 
rules (52). 

Dealing with dressed particles we consider only diagrams with thick fermion lines 
determined from the corresponding Dyson equations. Approximation (A.l) is how- 
ever based on the assumption that ImEf is much smaller than all other energies 
scales entering the problem. 

In particular in the limit ImE^ — in the fermion Green's functions one comes to 
the QPA, where the imaginary part of retarded Green's function becomes a delta 
function, e.g. for non-relativistic fermions 

ImGf ~ -^5[e + /^F-eJ--ReSP(e + /XF,p)] , Zp = 1/ {1 - dReT.^ /de){k.2) 

where the dispersion relation between e and p is implicitly given by 

ep-~eJ-+ReS|(ep.,p). (A.3) 
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Here the factor corrects for retardation effects, c.f. ref. [29], sucli that the sum 
rule (52) remains fulfilled. For simplicity we employ a quadratic p-dependence for 
ReS^ in terms of an effective fermion mass ep ~ /2m*p in the applications. 

Prom the dispersion equations one easily finds the corresponding relaxation times. 
E.g., for fermions supposing that ImE"^ is small and introducing T^^^y — l/25e, where 
we use that %l)p ~ exp(— i(ep- — i5e)t), one finds 

(A.4) 

This value would tend to infinity for ImSj^ — 0. In reality the fermion width de- 
termined by the value ImE^ is rather large even at sufficiently small temperatures. 
E.g. for nucleons, applying the QPA for the intermediate nucleon lines, it can be 
estimated as follows [54] 

ImEf ~ -\MfW1 - ep/epf + T^TrVe^], T <^ ep (A.5) 

with typically |MfP ~ at normal nuclear density [50]. This estimate shows that 
|ImS|^| comes into the order oiep already at sufficiently small temperatures T ~ |ei? 
and still increases for higher T. This defers the application of the QPA for a wide 
range of temperatures. 



B Renormalization of the Two— Fermion Interaction 

As an example we consider a theory where non-relativistic fermions interact via two- 
body potentials. A priori this theory has no bosons and the two-body interactions 
always connect two vertices of same sign, defining iV and iV'^'^ — {iV )\ 
while V~'^ — V'^~ — 0. Even if resumed to an effective four point interaction Qq~ 
according to eq. (63) much like Bruckner G-matrix, one has an effective interaction 
that connects only like sign vertices. 

Por the following we approximate Qo by a two-point function (as for instance in 
Permi-liquid theory, where the residual interaction is supposed to be local and ex- 
tracted from comparison with experimental data [20,10]), while Qq~^ = Qq = 0. 
We also suppose that Q^'^ and Qq~ interactions are particle-hole irreducible in the 
t-channel (vertical in (63), c.f. [20,10]). 

Starting from and ^o"*" one can completely bosonize the interaction in the 
standard way by resumming all intermediate particle-hole loop insertions 

A^^, = Gf,G',i, (B.l) 
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through the Dyson equation in two by two matrix form 



Q^Go + QooAeg. (B.2) 

In space-time homogeneous cases eq. (B.2) can be solved algebraically. First one 
defines a residual interactions through repeated A and insertions as 



1 - Go' A— ' 1 - g^^A++ 1 + Qo~A++ 

Using Qq^ = Qq^ = straight forward algebra yields 

g^- = zg^+A^-g;-,, g~+ = zg;-^A-^gxt, 

g++^ zgxt, g- = zg-:, (b.4) 

for the components of the full interaction Q, where 

z = (i - gXW-s;-sA-^y' (B.5) 

is a renormalization factor. 

The explicit dependence of g and g~^~ on Z can be moved to a renormalization 
of the loop 



g-+ = g-:A;+g+t, where A;+ = ZA-+. (B.6) 

This full interaction Q has bosonic features, just describing effective bosons, such as 
phonons, plasmons, sigma mesons, etc. Also the inclusion of real bosons, like pions 
in nuclear matter, is possible, giving rise to a picture, where pions couple to pionic 
particle- hole excitations, see [48,10]. These effective bosons can be taken on the same 
footing as all other effective quanta, the fermions or other bosons, see [43,39]. Thus, 
effective bosons also acquire a spectral function with width and the non-diagonal 
components of g~^ and g^~ are also Wigner functions. Consequently one comes to a 
theory of effective in-medium fermions interacting with effective in-medium bosons. 

Obviously the full (anti-)time-ordered g (^^^) depend on the in-matter den- 
sities G^^ and also via Z. However in an approximation where the value 
gQ^A^^g^ A^^ is small, one has Z ~ 1 for renormalization factor (B.5). So, one 
can simplify further and comes to a scheme like leading logarithmic approximation 
in quantum field theory, namely, a perturbation series over (or g'^~) neglecting 
corrections 1 + 0{g~^~A~'^) in each leading term. Such corrections are proportional 
to (or r^). Thus, one approximately has 

+ 

= ± ^ g;-sA^^g:r,t ^ J (b.t) 
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with like-sign effective interactions (63) (whereas in the general case one comes to 
(B.7) with renormalized A~f^ loops). 



C Contribution of More Complicated Diagrams 

Neglecting vertex corrections diagram (95) is given as 



-in-+i(g) = -iV^'VA^+iq) , where 



xiG~+{p)iG+-{p + k- q)ig-+{k) , 

with four-vectors q = {ojq.q), p = {e,k) and k = {u,k). Here is the " — h" 
interaction loop (B.7). For q <^ k ^ pp the integration in (C.l) over the p/c-angle 
givesP^ 



1 



Jo(e) = j dxC{p + k)G^'{p + k - q) 
-1 



,2m*^Ti -Uq + iT th((6 + ^)/(2T)) 
- (1 - ^e+.-.J^^ ^jr^, . 

Since | x | < 1 one has the following restrictions on p 

ep = p^/2m*p >eQ = {uj- k'^/2m*p - ujgfmp/2k^. (C.3) 

Further integrations can only be done in certain limits. For r,ujq <^T one can use 
Kadanoff-Baym ansatz (87), c.f. Appendix A, (A.l). In that case the e integration 
can be performed. One needs only the real part of expression (C.l) since the imag- 
inary part is cancelled by the corresponding diagram with opposite time ordering 
(opposite line sense). Thus, 



Ji = Re J deG++{p-q)G-^{p)Jo{e) 



■n a-n . 2m^vr^ u;^,-T\ 



where 



16 



using also q, k and p for \q\, etc. 
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^^th^^th ^ = 1 - 2n,_^^(l - n,+^) - 2n,+^{l - n,_^J (C.5) 



p — ^p—UJq ) 



For low fermion occupations u is about unity and we obtain with the help of relation 
(78) 



-U-+ = Ci(a;,) 1 ^^^n5<_JmEj5(u;,A;)J2(u;,A;)A;dA;da;, with (C.6) 



The integral in eq. (C.6) can be expressed through ^q"*" in QPA (eq. (82)), since 
pk^ m*T and one obtains 

xlm^^(a;, k)lniJ^{ijOq — k)k'^dkdu. 

This expression differs from the contribution of the corresponding QPA Feynman 
diagram calculated in ref. [45] only by the pre-factor Ci{ijjq) which is non-unit in 
our case of finite width P. In the QC hmit this expression (C.7) coincides with the 
n — 1 term in classical Langevin result. 

Diagram 

-\Aei (C.8) 



can be evaluated along similar lines considering the integrals for the left and right 
sub-loops 

J3= [{-l)iG—{p + k)iG-+{p)iG+-{p + k-q)'^--^{2^Y (C.9) 

J4 = /(-l)iG++(pi + ^-?)iG'-+(pi + ^)iG'+-(pi)^. (C.IO) 

Integration of (C.9), (C.IO) is done quite analogously to the previous cases. One may 
integrate over pfc-angle, then over e and ep in eq. (C.9) and over p^/c-angle, ei and 
epi in eq. (C.IO), respectively. One recovers the corresponding QPA form derived in 
ref. [45], however multiplied by the pre-factor Co(co'g), c.f. (95). 
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